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Abstract 


This  paper  presents  a  highly  desirable,  but  up  to  this  day  not  practical,  spherical  phased 
array  antenna  design  that  provides  hemispherical  coverage.  Consequently  it  can  be  used  for 
simultaneous  communication  and  control  with  multiple  satellites  and  air/space  surveillance  in 
the  entire  sky.  The  proposed  spherical  phased  array  antenna  structure  consists  of  a  number 
of  near-equilateral  triangular  planar  subarrays  arranged  in  an  icosahedral  geodesic  dome 
configuration.  This  novel  architecture  design,  ’’the  geodesic  dome  phased  array  antenna” 
preserves  all  the  advantages  of  spherical  phased  array  antennas  while  fabrication  is  based  on 
its  well  developed,  easily  manufactured,  and  affordable  planar  array  technology.  -  - 


f 


°This  work  was  supported  by  Air  Force  Research  Laboratory  under  Project  6.2  and  Aerospace  Corporation 
under  AFSPC/SMC  contract  FA8802-04-C-0001. 

B.  Tomasic  is  with  Air  Force  Research  Laboratory,  Sensors  Directorate,  Electromagnetics  Technology 
Division,  Antenna  Technology  Branch,  80  Scott  Road,  Hanscom  AFB,  MA  01731. 

J.  Turtle  is  with  Air  Force  Research  Laboratory,  Sensors  Directorate,  Electromagnetics  Technology  Divi¬ 
sion,  Antenna  Technology  Branch,  80  Scott  Road,  Hanscom  AFB,  MA  01731. 

S.  Liu  is  with  The  Aerospace  Corporation,  P.O.  Box  92957,  Los  Angeles,  CA  90009 

iii 


Contents 


Abstract  iii 

Contents  iv 

1  Introduction  1 

2  Spherical  Arrays  4 

2.1  Principle  of  Operation  4 

*  2.2  Why  Spherical  Arrays?  4 

2.3  Comparison  with  Other  Array  Architectures  5 

¥  2.4  Basic  Design  Parameters  8 

2.4.1  Active  Sector  8 

2.4.2  Element  Field  Pattern  9 

2.4.3  Element  Gain  10 

2.4.4  Array  Far-Field  10 

2.4.5  On-Off  Switch  10 

2.4.6  Excitation/Beamsteering  11 

2.4.7  Array  Gain  11 

2.4.8  Polarization  13 

2.4.9  Instantaneous  Bandwidth  17 

2.4.10  Mismatch  Losses  18 

2.5  Numerical  Simulation  19 

2.5.1  Computer  Code  19 

2.5.2  Calculated  Array  Patterns  20 

3  Geodesic  Dome  Phased  Array  Antenna  48 

3.1  Antenna  Architecture  48 

3.2  Geodesic  Dome  Structure  -  Background  ■  50 

9 

4  Geodesic  Dome  Antenna  Structure  *  62 

4.1  Global  Coordinate  System  62 

4.1.1  Counting  Icosa’s  62 

4.1.2  Counting  Sub-Icosa’s  62 

4.1.3  Class  I,  Method  1,  Icosahedral  Coordinates  62 

4.1.4  Class  I,  Method  1,  6v  Chord  Lengths  63 

4.1.5  Sub-Icosas  (Panel)  Types  ..  63 

4.2  Local  Coordinate  System  64 

4.2.1  Sub-Icosa  (panel)  in  Local  Coordinate  System  64 

4.2.2  Sub-Icosa  Base  67 

4.3  Reference  Coordinate  System  68 

4.3.1  Sub-Icosa  in  Reference  Coordinates  69 

4.3.2  Fitting  Equilateral  Triangle  Into  Sub-Icosa  Triangle  69 


iv 


4.3.3  Counting  the  Array  Elements  in  Reference  Coordinate  System  (x'NM,  y'NM)  72 

4.3.4  Location  of  the  Array  Elements  in  Reference  Coordinates  72 

4.3.5  Location  of  the  Array  Elements  in  Local  Coordinates  73 

4.3.6  Location  of  the  Array  Elements  in  Global  Coordinates  74 

5  Geodesic  Dome  Array  Beamsteering  and  Radiation  Pattern  116 

5.1  Element  Pattern  116 

5.1.1  Element  Pattern  in  Reference  Coordinates  116 

5.1.2  Element  Pattern  in  Local  Coordinates  116 

5.1.3  Element  Pattern  in  Global  Coordinates  117 

5.2  Far-Field  118 

5.2.1  On-Off  Switch  119 

5.2.2  Element  Voltage  Excitation  Coefficients  121 

v  5.3  Array  Gain  121 

6  Computer  Codes  128 

6.1  Code  1  128 

6.1.1  Input  128 

6.1.2  Output  128 

6.2  Code  2  129 

6.2.1  Input  129 

6.2.2  Output  130 

Appendix  A:  Element  Patterns  in  Source  (Local)  Coordinates  155 

A.l  Far-Field  of  Two  Crossed  Dipoles  155 

A. 2  Approximate  Expressions  for  Array  Element  Patterns  156 

Appendix  B:  Spherical  Array  Analysis  157 

t  B.l  Array  Model  157 

B. 2  Global  and  Local  Coordinate  Systems  157 

B.3  Coordinate  Transformation  Relations  ■  158 

*  B.4  Element  Patterns  in  Local  Coordinates  159 

B.5  Element  Patterns  in  Global  Coordinates  '  159 

B. 6  An  Example  160 

Appendix  C:  Maximum  Gain  Theorem  for  Conformal  Arrays  165 

C. l  The  Theorem  165 

C.2  Proof  166 

Bibliography  171 


v 


Chapter  1 
Introduction 


Within  the  next  decade  there  will  be  a  dramatic  increase  in  the  number  of  commercial 
and  military  satellite  constellations  providing  world- wide  telecommunication,  environmental, 
navigation,  and  surveillance  services  [1].  The  satellites,  in  LEO,  MEO,  and  GEO  orbits,  will 
need  a  number  of  ground  stations  to  serve  as  gateways  for  payload  message/data  routing 
and/or  as  control  centers  for  satellite  tracking,  telemetry,  and  command  (TT&C)  operations. 
A  typical  ground  station  will  require  one  or  more  large  high  performance  antennas  with 
hemispherical  coverage  for  communicating  with  LEO  and  MEO  satellites  during  the  brief 
interval  from  their  rise  above  the  horizon  until  they  set  below  the  horizon. 

Three  main  types  of  antennas  can  be  used  for  satellite  ground  stations:  mechanically 
steered  reflector  antennas,  electronically  steerable  arrays  (ESA),  and  microwave  lens  an¬ 
tennas.  Although  mechanically  steered  reflectors  are  most  widely  used  for  satellite  com¬ 
munications,  ESA  are  the  optimum  choice  for  this  application  because  they  can  provide 
multi-function,  multiple  simultaneous  beams,  have  high  reliability  with  graceful  degradation, 
require  low  maintenance,  and  have  low  life  cycle  cost.  With  rapidly  emerging  Microwave 
Monolithic  Integrated  Circuits  (MMIC)  and  multilayer  beamforming  and  beam  steering 
technologies,  large  ESA  are  viable  candidates  for  ground  station  antennas,  and  by  using 
commercial  off-the-shelf  components  they  could  be  also  affordable.  • 

Because  hemispherical  coverage  is  required  for  ground  station  antennas,  only  a  few  po¬ 
tential  array  architectures  are  suitable  [33].  These  can  be  classified  into  three  major  groups: 
(1)  multi-face  planar  arrays,  (2)  curved  surface  or  conformal  arrays  [12,  13,  14],  and  (3)  lens 
arrays/antennas.  In  [2,  3,  4]  we  present  design  trade-offs  for  these  array  architectures.  We 
have  shown  that  a  spherical  array  is  the  optimal  choice  for  ground-based  satellite  control 
antennas  in  terms  of  performance  and  cost.  It  is  a  natural  array  geometry  that  provides  scan 
independent  hemispherical  coverage  with  a  uniform  gain  pattern.  We  show  that  a  spherical 
array  requires  about  20%  fewer  elements  than  a  comparable  multi-planar  phased  array  [3]. 
It  is  also  shown  that  spherical  arrays  exhibit  significantly  larger  instantaneous  bandwidths 
and  have  lower  polarization  and  mismatch  losses  than  other  types  of  array  architectures. 
However,  fabrication  and  assembly  including  implementing  multi-beamforming  networks,  is 
much  more  difficult  for  curved  array  surfaces  than  for  planar  arrays  and  consequently  spher¬ 
ical  arrays  are  very  expensive.  This  is  the  main  reason  that  up  to  this  date  large  spherical 
phased  arrays  have  not  been  implemented  in  practice. 

A  novel  architecture  design,  ’’The  Geodesic  Dome  Phased  Array  Antenna”  that  preserves 
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all  the  advantages  of  spherical  phased  arrays  for  hemispherical  coverage  while  maintaining 
the  ease  of  fabrication  afforded  by  well-developed  planar  array  technology  was  invented 
and  proposed  in  [3]  by  Boris  Tomasic.  The  invention,  therefore  provides  desired  spherical 
phased  array  performance  using  affordable  planar  array  technology.  This  invention  is  also 
described  in  [3],  TT&C  Phased  Array  Antenna  Survey,  Attachment  3,  Section  4,  which 
was  sent  out  to  many  government  agencies  and  DoD  contractors,  and  cited  in  SBIR  Topic 
AF99-244.  The  geodesic  dome  array  was  also  analyzed  and  described  in  [4].  In  this  design, 
the  spherical  array  is  constructed  with  many  near-equitriangular  planar  subarrays  arranged 
in  an  icosahedra!  geodesic  dome  configuration.  This  "faceted”  dome  antenna  provides  a 
communication  function  for  SATOPS  and  a  radar  function  for  air/space  surveillance  with 
Ml  hemispherical  coverage  while  exhibiting  the  following  advantages  over  the  conventional 
pyramid-like  and  conformal  array  structures: 

1.  Keeping  the  architecture  locally  planar  and  globally  spherical  allows  the  array  to  be 
constructed  using  flat  (panel)  subarrays  which  are  assembled  to  form  the  hemispherical  dome 
structure.  This  is  technically  and  economically  viable  because  the  design  avoids  the  fabri¬ 
cation  complexity  associated  with  conformal  (curved)  arrays  since  the  subarray  fabrication 
is  based  on  well-developed,  easily  manufacturable  planar  array  technology.  Moreover,  the 
geodesic  dome  is  a  practically  realizable  structure  with  well-known  mechanical  design  and 
fabrication  techniques  [8,  9]; 

2.  The  Geodesic  Dome  antenna  preserves  all  the  advantages  of  spherical  phased  arrays 
[5,  6,  7]  such  as:  uniform  beams  over  a  hemisphere;  high  gain,  high  instantaneous  band¬ 
width,  low  mismatch  and  polarization  losses,  and  low  life  cycle  cost.  Above  all,  it  requires 
about  20%  fewer  radiating  elements  and  transmit/receive  (T/R)  modules  than  other  array 
configurations. 

We  have  designed  a  geodesic  dome  phased  array  antenna  that  meets  the  basic  require¬ 
ments  for  Space-Ground  Link  System  (SGLS)  and  Unified  S-band  (USB)  satellite  operation, 
[19]  to  [26].  The  antenna,  about  10  m  in  diameter,  will  support  two  independent  Tx  beams 
and  two  independent  Rx  beams.  It  consists  of  675  flat  panel/facets  (subdivision  frequency 
v  =  6).  One  of  the  panel  subarrays  was  fabricated  and  tested  by  AFRL/SNH.  Each  panel 
consists  of  a  multi-layer  beamforming  network  (BFN),  78  wide-band  circularly  polarized 
(CP)  radiating  elements  and  T/R  modules,  a  beamsteering  antenna  controller/computer 
(ACC),  and  software  that  will  integrate  the  geodesic  dome  antenna  with  control  computers 
in  the  existing  US  Air  Force  Satellite  Control  Network.  The  controller,  in  addition  to  per¬ 
forming  the  present  dish  antenna  functions,  will  be  capable  of  exploiting  all  ’’nice”  features 
of  a  phased  array  antenna  such  as  multiple,  simultaneous,  independent  beams,  fast  beam 
switching/steering,  adaptive  pattern  control,  etc.,  for  improved  SATOPS. 

The  report  is  organized  in  two  parts.  In  part  I,  we  present  an  extensive  theoretical 
analysis  of  spherical  arrays,  describe  the  geodesic  dome  structure  and  present  the  analysis 
for  excitation,  beamsteering  and  far-field  patterns  of  the  geodesic  dome  array  antenna.  In 
part  II  we  describe  the  system  aspects  of  the  Ml  dome  array  for  SGLS  and  USB  for  satellite 
TT&C,  and  radar  for  air  and  space  surveillance.  We  focus  on  design,  development  and 
demonstration  of  the  basic  building  block,  a  flat  panel  subarray,  for  SGLS  and  USB  TT&C 
satellite  operations  (SATOPS).  Specifically,  we  describe  design,  development  and  fabrication 
of  the  radiating  elements,  multi-layer  bemforming  network,  T/R  modules  and  beamsteering 
controller/computer.  Theoretical  and  measured  subarray  patterns  for  transmit  and  receive 
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beams  will  be  also  given. 

Part  I  of  this  report  consists  of  6  Chapters.  In  Chapter  2,  we  describe  conventional 
spherical  arrays  and  their  principle  of  operation.  Basic  design  parameters  such  as  element 
field  pattern,  array  far-field,  active  sector,  On-Off  switch,  excitation  and  beamsteering  are 
discussed  in  great  detail.  Also  we  derive  expressions  for  polarization  loss,  instantaneous 
bandwidth  and  mismatch  losses  in  spherical  arrays.  In  short,  this  chapter  treats  all  relevant 
design  aspects  of  the  spherical  arrays  in  multi-dimensional  trade  space,  that  will  serve  as  a 
basis  in  the  analysis  of  the  geodesic  dome  antenna.  In  Chapter  3,  we  give  background  on 
geodesic  dome  antenna.  The  geodesic  dome  structure  is  treated  in  Chapter  4.  The  array 
beamsteering  and  radiation  patterns  were  computed  in  Chapter  5,  while  in  Chapter  6  we 
describe  the  computer  codes  and  present  some  representative  numerical  data. 


Chapter  2 

Spherical  Arrays 


2.1  Principle  of  Operation 

The  general  configuration  for  spherical  arrays  consists  of  a  large  number  of  radiating  elements 
placed  on  a  spherical  surface,  as  schematically  shown  in  Figure  2.1,  The  basic  theoretical 
and  experimental  aspects  of  spherical  phased  arrays  are  given  in  [5,  6,  7],  For  a  given  beam 
direction  only  a  sector  of  the  array  elements  is  excited  while  all  other  elements  are  turned 
off.  The  active  sector  is  the  area  encompassed  by  a  cone  with  the  cone  angle  of  2amax 
degrees,  and  with  the  cone  axis  coinciding  with  the  antenna  beam  direction.  The  active 
part  (sector)  is  defined  as  that  part  of  the  sphere  that  is  turned  on  to  receive  (or  transmit) 
electromagnetic  energy.  Beam  scanning  is  accomplished  by  activating  different  sectors  of  the 
spherical  surface,  Figure  2.1.  All  array  elements  in  the  active  sector  are  phased  to  produce  an 
equi-phase  front  normal  to  the  desired  beam  axis.  Thus  each  array  element  must  have  an  RF 
On/Off  switch  and  a  phase  shifter.  The  radiating  elements  in  the  active  sector  are  assumed 
to  have  symmetry  in  the  plane  perpendicular  to  the  axis  of  the  cone  (antenna  beam),  and 
consequently  a  spherical  array  can  be  used  to  cover  the  hemisphere  with  practically  identical 
beams.  Thus  in  contrast  to  other  array  configurations,  which  suffer  from  beam  degradation  as 
the  beam  is  steered  over  wide  angular  regions,  spherical  arrays  can  provide  uniform  patterns 
and  gain  over  the  entire  sky.  In  order  to  provide  full  hemispherical  coverage  using  a  120 
degree  active  sector,  the  size  of  the  spherical  array  must  be  93%  of  the  Ml  sphere. 


2.2  Why  Spherical  Array? 

First  question  becomes  when  and  why  do  we  use  spherical  array.  To  answer  this  question  we 
present  a  simple  comparison  between  a  planar  and  spherical  arrays  in  terms  of  required  area, 
which  translates  into  number  of  array  elements  and  ultimately  the  array  cost.  We  begin  by 
considering  a  planar  array  of  circular  aperture  with  radius  a  and  area  Av  =  no2  as  shown  in 
Figure  2.2.  We  assume  that  the  array  is  required  to  scan  in  a  cone  with  angle  amax  measured 
from  z-axis.  The  peak  array  gain  at  8q  is 

Gp(e o)  =  47r^cos0o  (2-1) 

Similarly,  in  reference  to  Figure  2,3,  the  gain  of  the  spherical  array  is 
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CU  =  4,r^|pl  (2.2) 

where  the  projected  area  of  active  sector  is 

■^4  (o^max)  =  itR  sin  ocmax  (2-3) 

By  equating  (2.1)  and  (2.2),  the  radius  of  the  spherical  array  must  be 

R  =  — — — - yj  cos  ~9q.  (2.4) 

S1H  OLmax 

The  area  of  the  spherical  array  normalized  to  its  planar  counterpart  is  therefore 

AjOrnaJ  =  gj  ~  COS(^0  +  Omar)  CQg  gQ  (2.5) 

4p  Sin  Qtrmn-r 

Figure  2.4  plots  the  equation  (2.5)  for  amax  =  60°  and  90°  from  where  we  can  see  that  for 
scan  angles  up  to  roughly  75  degrees  the  planar  array  area  is  smaller  than  the  spherical,  and 
consequently  of  lower  cost.  For  scan  angles  greater  than  75  degrees  the  spherical  array  area 
becomes  smaller  than  the  planar  and  therefore  is  a  better  solution  for  this  application.  It  is 
worth  mentioning  that  in  this  simple  comparison  only  cos  6  scan  loss  is  taken  into  account.  In 
practice,  because  of  additional  mismatch  and  polarization  losses,  planar  arrays  are  not  used 
to  scan  beyond  60  to  70  degrees  from  broadside.  Thus,  spherical  arrays  complement  planar 
arrays  in  the  sense  that  for  small  scan  angles  (<  60°)  planar  arrays  are  more  efficient  than 
spherical  while  for  scan  angles  roughly  greater  than  60°  spherical  arrays  are  more  efficient 
than  planar  arrays. 

2.3  Comparison  with  Other  Array  Architectures 

As  shown  above  a  single  planar  array  cannot  be  used  for  wide  scan  angles.  In  this  case  we 
must  resort  to  the  alternative  array  configurations  such  as  conformal  arrays  or  geometries 
consisting  of  two  or  more  planar  faces.  Three  faces  are  the  minimum  required  for  hemisphere 
coverage.  Each  face  scans  throughout  a  portion  of  the  sector,  thereby  minimizing  the  max¬ 
imum  required  scan  angle  and  limiting  the  degradation  in  performance  caused  by  scanning 
the  beam.  In  general,  as  shown  in  [27]  and  below,  the  greater  the  number  of  faces,  the  better 
the  overall  antenna  performance.  In  [27,  28,  32,  34,  35]  it  is  also  shown  that  by  increasing 
the  number  of  faces,  less  elements  are  needed  for  a  given  array  performance.  (In  the  limiting 
case  a  multi  face  planar  array  geometry  converges  into  spherical  array  which  gives  the  best 
performance  with  minimum  number  of  array  elements.) 

In  the  remaining  part  of  this  section  we  show  that  in  the  limiting  case,  a  multi-face  planar 
array  geometry  converges  into  a  spherical  array  which  gives  the  best  performance  with  the 
minimum  number  of  array  elements.  For  convenience,  but  with  no  loss  of  generality,  we 
consider  a  simple  two-dimensional  (circular)  array  in  the  x-z  plane  with  nmax  number  of 
elements  as  shown  in  Figure  2.5.  We  also  assume  that  the  elements  are  uniformly  spaced  at 
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A/2  distance  and  are  equally  excited  and  phased  so  as  to  produce  a  main  beam  in  z-direction. 
The  element  gain  which  we  assume  to  be  the  same  for  all  elements  is 

g(9)  =  Gceii  cos  9  (2.6) 


where  6  is  the  field  angle  counted  from  the  local  element  normal.  In  (2.6),  the  unit  cell  area 
gain  is 


G  cell  — 


4tt  Aceu 
A2 


(2.7) 


with  Aceii  being  the  cell  area. 

Thus,  the  circular  array  peak  gain  is 

Tlrnax 

Gmax=  £dngn(9n-90)  (2.8) 

71—1 


where  the  Kronecker  delta  is 


i  if  |0„-flo|<f 
0  otherwise 


(2.9) 


and  0o  specifies  the  beam  direction  measured  from  the  2-axis.  The  array  gain  pattern  for  the 
selected  case  nmax  =  96  and  9q  —  0°  is  shown  in  Figure  2.6.  Figure  2.7  shows  the  peak  gain 
vs.  scan  angle  90.  As  expected,  because  the  array  geometry  is  independent  of  angle  6,  the 
gain  pattern  is  uniform  (independent  of  9)  and  the  peak  gain  is  also  practically  independent 
of  scan  angle  9q. 

A  polygon  approximation  to  a  circular  array  is  shown  in  Figure  2.8.  There  are  Nmax 
sides  with  N  =  1,2,...,  Nmax  being  the  sector  serial  number.  The  total  number  of  array 
elements  on  the  polygon  is  the  same  as  on  a  circle.  Now,  we  compute  the  peak  gain  of  the 
polygon-shaped  array.  The  peak  gain  is  no  longer  constant  vs.  scan  angle  9q,  but  periodically 
varies  between  Gmin  and  Gmax  to  be  determined  below. 

The  peak  gain  of  the  JVmQ2:-polygon  array  is 


7?  Nmax 

G(90)  =  Gceu- ~-2  £  SN(a)  cos  a  •  (2.10) 

™max  jv=l 

where  a  =  6^  —  8q  and  Kronecker  delta- is 

=  {  l  otherwise*  •  ^ 

We  distinguish  two  cases: 


(1)  Even-sided  polygon  with  two  sub-cases 

(la)  Number  of  sectors  given  by  Nmax  =  4 (m  +  1),  where  m  =  0,1,2,...  i.e.,  Nmax  = 
4, 8, 12, 16, 20, ....  In  this  case  the  maximum  peak  gain  is  at  90  =  ir/Nmax  and  minimum 
gain  occurs  at  9  =  0. 

(lb)  Number  of  sectors  given  by  Nmax  =  4m  +  6,  where  m  =  0, 1, 2, . . .  which  means 
Nmax  =  6, 10, 14, 18, 22, ....  In  this  case  the  maximum  and  the  minimum  gain  axe  at  90  =  0 
and  90  =  ■KfNmax ,  respectively. 


% 
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(2)  Odd-sided  polygon 

For  an  odd  sided  polygon  with  number  of  sides  Nmax  =  3  +  2m,  where  m  =  0, 1, 2, . . .  which 

is  Nmax  =  3, 5, 7, 9, _ ,  the  maximum  gain  is  at  0O  =  0°,  ir/Nmax  while  the  minimum  gain  is 

at  0O  =  7r/(2iVmoa.). 

In  all  cases  above  G{6q)  is  a  periodic  function  of  0o  with  period  A0jv  =  2tt /Nmax. 

Figure  2.9  shows  an  equilateral  triangular  array  where  only  one  side  (red)  is  active.  The 
gain  patterns  for  a  beam  pointing  at  0O  =  0°  and  0O  =  60°  are  shown  in  Figure  2.10  and 
Figure  2.11,  respectively.  Notice  that  at  0o  =  60°,  the  peak  gain  is  maximum  (about  20  dB) 
while  at  0O  =  0°  the  peak  gain  is  about  3  dB  lower.  This  is  explicitely  shown  in  Figure  2.12, 
which  shows  that  if  one  side  of  the  triangle  is  active  at  the  time  the  peak  gain  varies  vs.  scan 
angle  0o  by  about  3dB.  As  the  beam  scans  in  0.  from  0  to  360  deg  it  must  be  handed  over 
from  one  side  to  the  next  side.  For  single-sided  active  triangular  array  the  hand-over  angles 
are  0O  =  0°,  120°  and  240°. 

In  practice  a  large  variation  of  the  peak  gain  vs.  angle  is  undesirable  because  the  array 
size  must  be  determined  so  as  to  meet  minimum  gain  specs.  Thus,  in  the  triangle  case  it  is 
better  to  have  one  or  two  (depend  on  scan  angle)  active  sides  as  shown  in  Figure  2.13.  In  this 
case,  according  to  (2.10)  (case  (2)),  the  peak  gain  is  maximum  at  0O  =  0°,  60°  and  minimum 
at  0o  —  30°.  This  is  also  seen  in  Figure  2.14  which  shows  the  variation  of  the  peak  gain 
vs.  scan  angle  0O.  Notice  that  the  min  peak  gain  is  now  only  0.6  dB  below  the  maximum 
value.  The  array  gain  pattern  for  0o  =  0°  is  shown  in  Figure  2.15  which  is  a  substantial 
improvement  over  single-side  active  case. 

Figure  2.16  shows  a  square  array.  From  (2.10)  (case  (la))  the  minimum  peak  gain  is  at 
0O  =  0°  while  maximum  gain  is  at  0O  =  45°.  Figure  2.17  and  Figure  2.18  show  array  gain 
pattern  at  0O  =  0°  and  peak  gain  variation  vs.  0o- 

Figure  2.19  to  Figure  2.21  show  pentagon  arrays  and  respective  gain  patterns  for  0o  =  0° 
and  peak  gain  variation  vs.  0o-  In  this  case  the  maximum  peak  gain  is  at  0o  =  0°  and 
36°  while  minimum  peak  gain  occurs  at  0O  =  18°.  As  one  may  observe  from  Figure  2.14, 
Figure  2.18  and  Figure  2.21  as  the  number  of  polygon  sides  increase  the  peak  gain  variation 
vs.  0O  decreases  and  the  minimum  peak  gain  increases. 

Continuing  in  the  same  fashion,  Figure  2.22  to  Figure  2.24  show  the  hexagon  array 
with  corresponding  gain  pattern  for  0O  =  0°,  and  peak,  gain  variation  as  a  function  of  0o- 
Similarly,  Figure  2.25  to  Figure  2.27  show  the  same  for  a  12-sided  polygon  while  Figure  2.28 
to  Figure  2.30  show  the  same  for  a  24-sided  polygon.  And,  Figure  2.31  to  Figure  2.33  show 
the  same  for  a  48-sided  polygon. 

In  summary  we  observe  that  the  minimum  peak  gain  of  an  Nmai-polygon  array  increases 
with  Nmax  and  the  peak  gain  variation  decreases  as  a  function  of  scan  angle  90.  This  is  borne 
out  in  Figure  2.34  which  shows  the  minimum  and  maximum  peak  gain  variation  vs.  number 
of  polygon  sides,  Nmax.  It  is  seen  that  in  the  limiting  case  when  a  polygon  becomes  circle 
(Nmax  =  nmax),  the  maximum  and  minimum  peak  gain  curves  coincide  into  the  constant 
gain  line  of  a  circular  array.  At  the  same  time  the  minimum  peak  gain  of  the  polygon  array 
reaches  the  maximum  possible  value,  that  of  a  circular  array. 

Thus,  we  conclude  that  if  omni  scan  coverage  is  required,  the  circular  array  is  the  optimal 
geometry,  since  it  offers  the  best  performance  in  comparison  to  polygon-type  arrays  with  the 
same  number  of  radiating  elements. 


In  summary,  the  spherical  array  offers  the  optimal  geometry  for  hemispherical  coverage. 
Beside  providing  a  uniform  beam  in  a  hemisphere  it  also  requires  minimum  number  of  radi¬ 
ating  elements  in  comparison  to  alternative  array  architectures  consisting  of  multiple  planar 
arrays  or  other  conformal  arrays. 

Below  we  present  a  simple  analysis  and  compare  a  spherical  array  with  a  conventional 
three-face  pyramid  structure.  Figure  2.35  shows  the  array  where  each  face  consist  of  circular 
planar  arrays  of  radius  a.  To  cover  a  hemisphere  each  individual  planar  array  must  scan 
about  64  degrees  from  its  local  normal  (broadside). 

If  area  of  the  single  face  is 

Ap  =  7r  a2  (2.12) 

then  the  total  array  area  is 

Aa  =  3Ap.  (2.13) 

Thus,  if  we  set  in  (2.5)  Ap  =  AAj 3  and  for  sake  of  simplicity  6$  =  7r/3  we  get 


A  [ctmax)  —  0AA 


2  .  1  C0S(g  +  OCjnax)  7T 


sin  a„ 


cos  — . 


(2.14) 


We  see  that  if  amax  =  90°  then  A(amax)  =  0.62.4a  and  if  amax  =  60°  then  A(amax)  =  2/3Aa. 
Thus,  a  spherical  array  requires  about  33%  less  array  elements  than  the  three-face  pyramid. 

Another  candidate  array  configuration  for  hemispherical  coverage  is  the  Dome  antenna 
[29,  30].  The  pros  and  cons  of  the  dome  antenna  are  discussed  in  [31]. 


2.4  Basic  Design  Parameters 

2.4.1  Active  Sector 

Once  we  have  decided  on  a  spherical  array,  the  question  becomes  how  big  is  the  active  sector 
or  what  is  the  optimum  value  of  amax.  The  array  gain  is 


G (ftmta)  —  {kR)  SU1  Ot-max 

(2.15) 

with  the  beamwidth 

BW  =  58°——^ - ■. 

kRsmamax 

(2.16) 

For  the  fixed  G/BW,  i.e. , 

kR  sin  amax  =  58  =  C 

BW 

(2.17) 

we  see  that  the  active  area  must  be 

A(amax)  =  ^ C2f(amax ) 

(2.18) 

where 

r/  \  _  1  Q-max 

““  •  o 

(2.19) 

Sill  Qmax 

As  seen  in  Figure  2.37  the  required  active  array  area  is  the  smallest  when  amax 

=  90°.  For 

practical  reasons  amax  is  usually  between  60  and  70  degrees. 


\ 
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2.4.2  Element  Field  Pattern 


The  fax  electric  field  produced  by  a  single  element  is  originally  expressed  in  a  local  coordinate 
system  (rnm,  9nm,  <f>nm)  with  the  origin  at  element  ( nm )  as  follows 


®nm(rnmi  ^nm) 


g j(ut-krnm) 


(2.20) 


where 


®nm(^nmi  <Anm)  —  ^nm^6nm{9nmt  $nm)  "h  finin')  . 


(221) 


where  egnm  and  e^m  are  the  complex  element  pattern  factors  of  the  two  field  components. 
Basically  they  represent  the  field  produced  by  a  single  element  in  a  match-terminated  spher¬ 
ical  array  environment  where  6  is  the  phase  difference  between  the  two  components  of  the 
field.  In  general,  the  field  is  elliptically  polarized.  For  a  circularly  polarized  field  egnm  =  e^,nTO 
and  S  =  ±7r/2  where  the  plus  sign  corresponds  to  RHCP  and  the  minus  sign  to  LHCP  of 
the  outgoing  wave. 

We  consider  an  elementary  case,  the  rough  approximation  for  the  pattern  of  two  CP 
crossed  infinitesimal  dipoles,  which  from  Appendix  A.4  is 


where 


$nm)  —  Eq  \J COS  6nm  6  ^nm  i  j<Pnm) 


(2.22) 


(2.23) 


This  pattern  can  be  expressed  in  global  coordinates  (r,  9 ,  fi)  using  coordinate  transformation 
relation  (B.12)  as  described  in  Appendix  B.6,  i.e., 


where 


enm(&,<t>)  =  —Eq  a/cOS ~dZl-j=(6  ±  j0)e±>nm(0,0)eiVnme±^nm 


cos  9nm  =  sin  an  sin  6  cos(0  —  (5nm) .+  cos  an  cos  6 

cos  an  sin  9  cos  (fi  —  (3nm)  —  sin  an  cqs  9 

cot  (prim  —  ;  21  •  7~1  o  \ 

sin  9  sm(fi  -  Pnm) 


ipnm  =  kR  cos  9nm  =  kR[ sin  an  sin  9  cos(fi  -  /3nm)  +  cos  an  cos  9] 


(2.24) 


(2.25) 

(2.26) 


(2.27) 


tan  vnm(9,fi)  = 


sin  an  sin(fi  ftnni) 

cos  9  sin  an  cos(<fi  —  f3nm)  —  sin  9  cos  an ' 


(2.28) 
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2.4.3  Element  Gain 

Element  gain  is  defined  as 


_  //,  jl  \  _  enm(r)  •  e’m(r) 

ynm\unmi  /*  p 

SO-* Dnm 

(2.29) 

where 

p  _  |  Cnm  | 

4wL 

and  Co  =  376.7  0  is  the  free  space  impedance.  Thus, 

(2.30) 

4^r  2 

9nm{$nm<}  $nrn)  ~  ~7~  |  $nrn}  |  * 

Co 

(2.31) 

Specifically,  if  enm 

above  Is  that  of  (2.24),  we  get 

9nm  ~  GqzH  COS  0nm 

(2.32) 

where 

rt  4 irAcdl 

Ucell  =  ^2 

(2.33) 

is  the  unit  cell  area  gain. 


2.4.4  Array  Far-Field 

The  expression  for  the  radiation  pattern  produced  by  the  spherical  array  is  obtained  after 
combining  the  fields  produced  by  the  individual  antenna  elements  at  the  far  field  point. 
Thus, 

ej(uit-kr ) 

E(r;t)  =  E(M) -  (2-34) 

r 

where 

E(0,  <f>)  =  9E6(9,  <l>)  +  <^(0,  $  =  £Enm(0,  <f>)  (2.35) 

„  nm 

and 

E nm(0,  <f>)  =  Snm  cnrn(0Q,  60)  enm(6,  <f>).  (2.36) 

Here,  5nm  are  on-off  switches,  enm  are  complex  element  patterns,  and  cnm  are  element  voltage 
excitation  coefficients. 

2.4.5  On-Off  Switch 

As  already  mentioned,  6nm  in  (2.36),  is  an  On-Off  switch.  This  switch  determines  which 
element  will  be  on  or  off  for  a  particular  beam  direction.  For  example,  if  we  define  an  active 
array  region  within  cone  with  axis  in  (9q,  do)  beam  direction  and  with  the  cone  angle  amax, 
then 

r  _  J  1  if  7nm  —  Otmax  fry 

dnm  ~  \  0  otherwise  {  ) 


where  7nm  is  the  angle  between  the  axis  of  the  cone  and  the  radial  line  passing  thru  the 
element  ( nm ).  Since 

cos  7nm  =  f  •  pnm  (2.38) 

where  pnm  is  the  radial  vector  specifying  the  location  of  the  (n,  m)th  element,  see  Figure  5.3, 
we  can  write 

cos  7„m  =  sin  d0  cos  </>0  sin  an  cos  j3nm 
+  sin  do  sin  sin  an  sin  (3nm 

+  cos  $o  cos  Oin.  (2.39) 

2.4.6  Excitation/Beamsteering 

The  excitation  coefficients  axe 

Cnm{0o,<po)  =1  CnmiOo^o)  |  (2.40) 

where  |  c„m  |  depends  on  the  desired  aperture  illumination  function.  The  array  elements  are 
phased  so  that  the  pattern  produced  by  the  array  has  a  maximum  along  the  (0q,  <po)  direc¬ 
tion.  This  means  that  Xnm(d o,  ^o)  must  be  of  the  opposite  sign  of  the  phase  of  enm(60,  (f>o)- 
Specifically,  for  enm(#o,  4>o)  of  (2.24),  the  phase  of  the  excitation  coefficients  is 

Xnm(.d o,  )  =  ~\}Pnm{d 0,  <^o)  i  Vnm[6oi  00  )  i  0nm(^ 0>  0o)]  (2.41) 

chosen  so  that  Re[E<jmrn{6o,  <f>0)  exp(jut )]  =  0  when  uit  =  0.  Thus  when  ut  =  0  all  the 

element’s  instantaneous  fields  point  in  the  same  (0)  direction. 


2.4.7  Array  Gain 

The  array  gain  is 


r,„  E(r)-E-(r) 

G(M)~  („Fo  1 

(2.42) 

where 

p  _  1  Cnm  | 

°  47rr2 

(2.43) 

Thus, 

( q  j.\  _  4tt|  ^2nm  CnmGnm(B ,  (j)')  | 

}  0)  /•  |  |2 

SO  2^mn  I  Cnm  j 

(2.44) 

and  consequently,  the  peak  gain  is 

r<{ a  A  \  ^  Enm  1  cnm  |  |  enm(^0j  00 )  |] 

^{do,(po)=  r  .  |2 

(>0  2-mm  1  Cnm  | 

(2.45) 

For  an 

element  pattern  as  given  by  (2.22)  we  get 

s~i(  r\  i  \  _  /-t  [5Znm  I  Cnm  ||  ^/COS  9nm  |] 

00  )  celt  1  12 

2^i  nm  1  ^nm  | 

(2.46) 
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The  expression  (2.45)  is  maximum  if  we  choose  the  excitation  coefficients  Cnm  according 
to  the  Maximum  Gain  Theorem,  (C.2),  which  is  for  convenience  presented  in  Appendix 
C.  Namely,  in  conformal  arrays,  in  contrast  to  planar  arrays,  maximum  gain  is  achieved 
not  when  elements  are  uniformly  excited  but  when  the  magnitude  of  the  voltage  excitation 
coefficients  is  proportional  to  the  magnitude  of  the  complex  element  pattern  (field)  values 
in  the  rnain  beam  direction  (6o,<j>o)i  and  the  phase  of  the  voltage  excitation  coefficients  is 
such  as  to  compensate  for  the  different  propagation  paths  and  element  phase  pattern  values 
in  {Oo,  (bo),  and  the  orientation  of  the  radiating  elements.  In  short,  (2.45)  is  maximum  when 

|  Cnm  |0C|  6nm  (*>,*>)!  (2-47) 

while  the  phase  of  Cnm  must  be  that  of  (2.41).  Thus,  for  a  spherical  array  of  radiators  with 
an  element  pattern  of  (2.22),  the  maximum  gain  can  be  achieved  if 

I  Cnm  1=  \] cos  6nm.  (2.48) 

As  an  illustrative  example,  below  we  evaluate  the  array  gain  for  (1)  "equal-amplitude” 
excitation,  i.e.,  |  c„m  |=  1  and  (2)  "maximum  gain  excitation”  as  given  by  (2.48).  Because 
the  radiation  pattern  in  a  spherical  array  is  independent  of  the  beam  pointing  direction,  with 
no  loss  of  generality,  we  assume  that  all  the  elements  in  the  active  sector  are  phased  so  that 
the  peak  of  the  beam  points  in  z— direction,  i.e.,  do  —  <f>o  =  0. 

In  case  (1),  from  (2.46) 

Gu(0,  0)  =  Gcell  "atE  \! C0S  ®nm\2  (2.49) 

iva  nm 

where  Na  is  the  number  of  elements  in  the  active  array  sector,  while  in  case  (2) 

Gmax(0,  0)  =  GcM  £cOS0nm  =  £<7nm(0,0).  (2.50) 

nm  nm 

The  above  sums  can  be  evaluated  by  replacing  Dn  Dm  with 

/*2tT  fOLmaa;  _  ,  . 

/  /  -  -  R2  sin  a  da  d/3.  (2.51) 

Jo  Jo 

Similarly,  in  (2.49),  Na  is  replaced  by  the  active  area  i.e., 

Na  -*  2itR2(1  -  cos  amax).  (2.52) 


For  amax  =  7r/3  the  ratio  between  (2.49)  and  (2.50)  is 


<?max(0;  0) 

G„(0,0) 


1.0086 


(2.53) 


Thus,  in  case  of  "maximum  gain  excitation”  the  gain  is  about  0.04  dB  higher  than  if  the 
array  elements  were  excited  uniformly.  This  increase  in  the  gain  is  insignificant  for  any 
practical  application. 


* 
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2.4.8  Polarization 

The  far-field  pattern  of  a  single  element  in  array  is 


E nm{9, 4>)  =  Cnm{6 o,  <po)  enm(0, 4>)  (2-54) 

For  the  sake  of  simplicity  in  the  previous  analysis  we  assumed  that  the  array  element, 
radiates  a  CP  wave  in  all  directions,  as  expressed  by  (2.22).  For  a  more  realistic  (physical) 
element,  (2.22)  is  modified  according  to 

enm(0nm,  0nm)  =  Eq  ^/cOS^Z^  e±j*nm -j=(dnm  COS  Qnm  ±  J0nm)  (2.55) 

which  resembles  the  far-field  radiated  by  two  CP  infinitesimal  dipoles  A/4  above  an  infinite 
ground,  Appendix  A.  The  field  is  elliptically  polarized  everywhere  except  at  the  local  zenith 
{6nm  =  0)  where  it  is  CP. 

We  express  this  field,  see  Appendix  B.6  for  details,  in  global  coordinates  as 


enm(M)  =  -[0  I  I  ±jj>  |  e*nm(M)  I 

|  0)  |=  E0  \J COS  0nm  \  CLnm  COS  0nm  i 

|  &(finm  ($  ?  0)  |=  fK  E0  COS  @nm  |  & nm  ^  COS  Q> 


'nm.  e±o4>  nm 


V2' 

tani?nm  = 
tan  (pnm  — 


®nm  COS  &nm 
&nm  COS  OjiTfi 
Qnm 


(2.56) 

(2.57) 

(2.58) 

(2.59) 

(2.60) 


where  anm  and  bnm  are  elements  of  the  rotation  matrix  given  by  (B.9)  and  (B.10). 
We  rewrite  (2.56)  in  the  form 


=  -[»  |  eenm(eA)  I  ± &  |  (2.61) 

where  for  convenience  we  set  ipnm  —  unm.  In  (2.54),  the  voltage  excitation  coefficients  are 

cnm(0o,  0o)  =|  Cnm{60)  0o)  I  (2.62) 

where,  as  already  mentioned,  to  point  the  beam  in  the  (0O, 0o)  direction,  Xnm  must  be 

Xnm{0 0,  00  )  =  ~bPnm(0 0,  00 )  ±  Vnm{6{ 0,  00 )  ±  0nm(#O,  0o)]-  (2-63) 

To  observe  the  polarization  effects  we  find  the  instantaneous  far  field  of  a  single  (nm) 
element  in  the  (0q,  0o)  direction: 


where  according  to  (2.54)  using  (2.62)  and  (2.63) 

Enm(0o,0o)  =  -  |  CnmiQoM  |  [0  |  e9nm(&o,  <h)  1  I  e^mOMo)  |] 

(2.65) 

Substituting  (2.65)  into  (2,64)  yields 

Enm($0s  ^Oi  t)  =  |  Cnmifi (h4*o)  I  {@  |  &8nm.{$Ch  00  )  j  COs[(l?nm(00,  0o)  ''Pnmi^Qh  00  ))  i  cut] 

=F  0>  I  e*nm(0o,0o)  [sinwf} 

(2.66) 

Thus  in  the  plane  perpendicular  to  the  array  beam  direction,  the  instantaneous  electric  field 
Enm  traces  an  ellipse  with  tilt  angle 


1  sin(2^nm)  COS  ig  * F  (^nm  fPnm\ 

(2.67) 

.  2  cos(27nm) 

where 

,  —1  &4>nm  | 

7nm  =  tan  1  p  . 

1  &Qnm  | 

(2.68) 

The  axial  ratio  is 

Rnm  ”  COt  €nm 

(2.69) 

where 

^nm 

X  7T 

=  -  sin_1{sin(27nm)  sin[±-  T  (4  -  'Am)]}- 

(2.70) 

For  simplicity  we  assume  that  we  point  a  beam  in  2-direction,  (0Oi0o)  =  (0,0),  with  array 
active  elements  located  within  a  cone  angle  amax  =  60°.  Then  only  the  top  element,  (nm)  = 
(00)  with  ao  =  0,  0oo  =  0  designated  as  the  reference  element,  will  be  CP  while  all  others 

will  be  elliptically  polarized.  Figure  2.39  shows  in  red  the  polarization  circle  of  the  reference 
element  in  the  z=constant  plane.  As  an  illustrative  example,  the  Eoq  vectors  for  ut  = 
0,  7r/2,  7r  and  ut  =  2>n/2  are  shown.  In  the  same  figure  we  also  show.the  polarization  ellipse 
of  the  (n,  m)  =  (38, 0)  element  with  coordinates  (an  =  30.4°,  @no  =  0°);  thus  the  element  is 
located  in  the  x-z  plane.  The  instantaneous  electric  field-  vector  Enm  is  shown  in  blue  for  the 
same  ut  parameter  values.  Notice  that  for  ut  =  0, 7r/2,  it,  and  3tt/2  both  vectors  are  in  sync 
while  for  other  ut  values  they  point  in  slightly  different  direction  resulting  in  a  reduction  of 
the  array  gain,  or  a  so  called  polarization  loss.  In  regions  0  <  ut  <  n/2  and  ir_<ut<  3ir/2 
the  Enmis  leading  E0q  while  in  regions  7r/2  <  ut  <  n  and  3ir/2  <  ut  <  2ir  Enm  is  lagging 
E00 .  The  element  phasing  was  chosen  according  to  (2.63)  so  that  at  ut  =  0  all  element 
instantaneous  vectors  align  with  the  top  element  (0,0),  as  can  be  seen  in  Figure  2.39. 

The  axial  ratio  and  orientation  of  the  ellipse  are  indicative  of  the  element  pattern  relation 
(2.55)  and  the  element  position  in  the  array.  Below  we  elaborate  on  the  ellipse  characteristics: 

1.  The  radius  of  the  red  circle  in  Figure  2.39  is  E0/-\/2,  see  (2.55)  where 


fbrAceti 


E0  =  ^t  =  ^ 


(2.71) 


i 

I 


i 

l 
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We  assume  a  typical  value  for  the  cell  area  Aceu  —  A/2  x  A/2  which  gives  an  area  cell  gain 
Gceii  =  tt.  Thus  in  (2.71),  E0  =  1.2533. 

2.  The  ellipse  axis  depends  on  element  location.  For  (n,  m)  =  (38, 0)  element,  the  major 
axis  of  the  polarization  ellipse  is  smaller  than  the  diameter  of  the  red  circle,  and  is  due  to 
the  y/cos  6nm  factor  in  (2.55). 

3.  The  ellipse  axial  ratio  is  determined  by  the  factor  (6nm  cos  9nm  ±  j4>nm)  in  (2.55). 
Thus  the  axial  ratio  is 

R  = - v-  (2.72) 


cos  8n 


or  in  our  case,  i.e.,  the  ^-directed  beam 


R  = 


cos  an 


(2.73) 


For  element  (38, 0)  shown  in  Figure  2.39,  R  =  1/  cos  30.4°  =  1.15. 

4.  The  ellipse  tilt  angle  can  be  computed  from  (2.67)  and  for  all  elements  lying  in  (+x)-z 
plane  it  is  90°. 

5.  Note  also  that  at  ut  =  0  both  vectors  point  in  the  negative  8  direction.  Below  we  give 
a  short  explanation:  It  is  assumed  that  at  ut  =  0  the  single  element  has  only  Egnm  in  the 
local  coordinate  system.  Because  the  m  =  0  element  is  located  in  the  x-z  plane,  ESnm  lies 
also  in  the  same  plane  as  shown  in  Figure  2.43.  The  elements  are  phased  so  as  to  produce  a 
beam  in  the  z-direction  (0o,0o)  =  (0,0),  and  therefore  the  corresponding  local  observation 
angles  are  (0Onm,  <Ponm)  =  (an,  — tt).  Transforming  this  field  into  global  coordinates  gives  -Eg, 
as  illustrated  in  Figure  2.43. 

Similarly,  Figure  2.40  shows  the  polarization  ellipse  for  the  (n,  m)  =  (76, 0)  element.  This 
element  is  located  at  the  edge  of  the  active  sector,  more  specifically  at  (a„,  (3nm)  =  (60.7°,  0°). 
Comparing  this  figure  with  Figure  2.39,  we  notice  that  the  ellipse  is  much  smaller  because 
of  the  .y/cos  an  factor  and  the  axial  ratio  is  about  MB.  This  is  not  great  but  still  acceptable 
performance. 

As  m  increases  the  polarization  ellipse  starts  rotating  in  the  /3nm  direction  and  makes  a 
full  circle  as  Pnm  sweeps  from  zero  to  27T.  A  snap  shot  in  time  of  Enm'  for  the  element  (nm)  = 
(38, 16)  with  an  =  30.4 °,/3nm  =  29.2°  and  the  corresponding  ellipse  for  ut  —  0,  7t/2,  i r,  3tt/2 
is  shown  Figure  2.41.  Because  of  different  ellipse  tilt  angles  as  a  function  of  element  position 
(o;n,  0nm)  the  E nm  vectors  do  not  rotate  in  sync  vs.  ut.  In  addition,  Enm  at  some  reference 
time,  for  example  ut  =  0,  does  not  coincide.  This,  however,  can  be  compensated  for  by 
imposing  the  proper  phase  of  the  excitation  coefficients  according  to  (2.63).  In  this  case,  at 
the  peak  of  the  beam,  Enm  will  align  with  E0o  at  four  distinct  points  on  the  ellipse:  utx  =  0, 


Utn  =  7T, 


-1  iflnm!  CQS^nm 

¥  sin(i?nm  -  ipnm) 


utz  —  tan 


(2.74) 


and  uti  =  utz  +  7r-  For  other  ut' s,  the  instantaneous  vector  fields  Enm  from  different 
array  elements  won’t  be  exactly  in  sync,  and  therefore  there  will  gain  degradation  due  to 
polarization  effects.  For  this  particular  element  utz  =  35.0474°  and  ut±  =  215.0474°,  see 
(2.74).  Also  notice  that  in  the  regions  0  <  ut  <  utz  and  7 t  <ut  <  ut$,  638,16  is  leading  Boo 
while  in  regions  utz  <ut<ir  and  2ir  <  ut  <  ut4  the  E38,i6  is  lagging  Eqq. 
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Figure  2.42  shows  the  same  as  Figure  2.41  except  here  the  element  is  (n,m)  =  (76,29). 
This  element  has  the  following  coordinates:  an  =  60.7°  and  /3nm  =  30.7°.  Furthermore, 
the  four  points  in  time  where,  the  instantaneous  electric  field  vectors  align  with  the  top  CP 
element  (0, 0)  are  two  fixed  values  ut  —  0,  tt  and  two  element  dependent  values  utz  —  43.1548° 
and  cut4  —  223.1548°. 

Below  we  evaluate  the  polarization  efficiency  (loss)  in  spherical  arrays.  We  define  the 
polarization  efficiency  in  spherical  arrays  as  the  ratio  of  the  peak  gain  with  elliptically  polar¬ 
ized  elements  to  the  peak  gain  of  the  array  when  all  elements  are  circularly  polarized.  Thus, 
using  (2.45)  for  the  array  gain,  we  get 


p=m 


Cnm 


»(%  $o)  I]2 


[E* 


I  II  «£(*.*»)  IP 


(2.75) 


where  enm  is  given  by  (2.55)  and  is  given  by  (2.22).  Under  maximum  gain  condition* 
i-e.,  |  Cnm  |=|  e(60,<p0)nm  |  and  |  |=|  e(0o,  <pQ)°*  |  and  using  transformation  (2.51), 

(2.75)  can  be  written  in  the  form 


which  gives 


For  a„ 


[fomax  sin  a  cos  a  ^  da}2 

P  [70“m“  sin  a  cos  a  da]2 

(2.76) 

2  [22/3  —  (2  —  sin2  amax)3^2}2 

P  9  sin4  amax 

(2.77) 

0.809  and  consequently  the  polarization  loss  is 

LP  —  — lOlogp  =  0.92  dB. 

(2.78) 

Evaluating  the  polarization  loss  in  the  planar  array  case,  (2.75)  is  still  valid  and  because 
all  elements  point  in  the  same  direction  it  can  be  further  reduced  to 


P- 


enm(#0,  4>o)  |2  [En  m  |  Cnm 


a£SM>.*>)lJE»,l«l  P' 

In  planar  arrays  maximum  gain  is  obtained  when  |  cnm  |=  1.  Therefore, 

_  |  ®nm(^Oi^o)  P 
|2' 

Substituting  (2.55)  and  (2.22)  into  (2.80)  yields 

1  +  cos2  60 
P= - o - 


(2.79) 


(2.80) 


(2.81) 


which  at  60°  off  broadside  scan  is  p  =  0.625,  and  polarization  loss  Lp  =  2.04 dB.  Thus,  the 
polarization  losses  of  the  spherical  array  are  about  half  of  the  equivalent  planar  array  at  60° 
off  broadside  scan. 
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2.4.9  Instantaneous  Bandwidth 

In  planar  arrays  steered  by  phase  shift  rather  than  by  time  delay  the  beam  will  scan  (squint) 
as  the  frequency  changes.  This  beam  squint  with  frequency  results  in  a  decrease  in  gain 
because  the  beam  is  no  longer  pointing  in  the  desired  direction.  Expressions  for  the  signal 
bandwidth  of  planar  phased  arrays  are  derived  in  [18,  12]  where  the  bandwidth  criteria  is 
based  on  unacceptable  gain  loss  due  to  change  in  scan  angle  with  frequency.  For  example, 
if  we  define  the  frequency  limits  in  which  the  beam  squints  by  1/4  of  the  beamwidth,  and 
looses  about  3/4  dB  gain,  the  resulting  bandwidth  is  given  by 


A  /  _  A 
/  2  L  sin  $o 


(2.82) 


For  0o  =  60°,  this  expression  reduces  to 


Bandwidth  (%)  «  Broadside  Beamwidth  (degrees).  (2.83) 


Thus  beam  squint  with  frequency  is  the  fundamental  mechanism  that  limits  the  band¬ 
width  of  planar  phased  array.  In  spherical  arrays  there  is  no  beam  squint  since  the  beam 
steering  principle  is  different  than  in  the  planar  arrays.  However,  in  spherical  arrays  there 
is  a  beam  broadening  effect  which  causes  gain  loss  as  the  frequency  changes.  Consequently, 
the  bandwidth  criteria  for  planar  arrays  developed  in  [18]  do  not  apply  in  a  spherical  array 
case. 

Below  we  develop  the  bandwidth  criteria  for  a  spherical  array  aperture.  The  criteria  is 
based  on  the  requirement  that  all  elements  in  an  active  array  sector  contribute  constructively 
to  the  total  far  field  in  the  beam  direction.  This  means  that  the  difference  in  phase  between 
the  center  element  and  the  element  at  the  edge  of  the  active  sector  is  not  greater  than  7r/2 
as  the  frequency  changes  from  /0  to  /.  Thus,  with  the  help  of  Figure  2.38  we  write 

j  Aj()-ft(l  COS  Qtfnax')  kR(\  COS  Otjnax  \  •  (2.84) 


From  here 


fo  f  |  (1  cos  Otmax)  —  ^  • 


(2.85) 


Setting  A /  =  2  |  /o  —  /  |  and  since  c  =  2irfo/ko,  the  instantaneous  bandwidth  of  a  spherical 


array  is 


A  /  _ _ 7T _ _ 

fo  k0  ft  ( 1  COS  OLjmxx ) 


(2.86) 


This  criteria  gives  about  1  dB  loss  at  the  edge  of  the  band.  As  an  illustrative  example,  we 
consider  a  spherical  array  of  5  meter  radius  with  an  active  sector  cone  angle  amax  =  60°  and 
/  =  2  GHz.  In  this  case  from  (2.86),  A///0  =  3%  or  A/  =  60  MHz.  Comparing  this  with  an 
equivalent  planar  array  with  the  same  projected  circular  aperture  of  radius  ft  sin  amax,  and 
with  a  broadside  beamwidth  of  1°;  according  to  (2.83),  the  signal  band-width  of  the  planar 
array  is  approximately  1%.  In  this  case  at  60°  off-broadside  scan  a  spherical  array  has  about 
three  times  greater  signal  bandwidth  than  an  equivalent  planar  array  . 

Figure  2.49  and  Figure  2.50  show  a  <p  =  0  cut  radiation  pattern  of  a  spherical  array  and 
an  equivalent  planar  array  with  a  beam  pointing  in  (do  =  60°,  0O  =  0°)  direction  for  two 
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frequencies  f0  =  2  GHz,  and  fo  4- 10  MHz  corresponding  to  1%  signal  bandwidth.  Prom 
Figure  2.50  notice  that  at  /0  the  planar  array  has  3  dB  scan  loss  while  the  spherical  array 
has  a  scan  independent  gain  pattern.  As  the  frequency  changes  from  2.000  GHz  to  2.010 
GHz  a  planar  array  exhibits  a  beam  shift,  and  an  additional  1  dB  loss,  while  a  spherical 
array  keeps  pointing  its  beam  in  the  same  direction  with  only  about  0.1  dB  loss  in  gain. 


2.4.10  Mismatch  Losses 

Let  anm  and  bnm  be  the  voltages  of  the  incident  and  reflected  signals  at  the  array  element 
terminals.  The  total  voltages,  also  called  array  excitation  coefficients,  c„m  are  then 

Cnm  ~  Gum  4"  bnm  (2.87) 

In  phased  arrays  the  incident  and  reflected  voltages  are  related  via  scan  dependent  active 
reflection  coefficients, 

bnm  =  Lnm  (^0i  00 )  ®nm  (2.88) 

where 

Tnm(0o,  <Po)  =|  Tnm(e0, 4>o)  I  '  (2.89) 

In  spherical  arrays,  as  in  planar  arrays,  the  excitation  coefficients  c„m  directly  affect  the 
array  gain  pattern,  and  are  usually  determined  according  to  the  desired  pattern  whether  in 
regard  to  the  sidelobe  level  or  maximum  gain  condition.  With  known  cnm’s,  the  incident 
signal  voltages  must  be 


Mismatch  losses  are  defined  as  the  ratio  of  the  input  or  available  power  to  the  power  radiated 
by  the  antenna,  i.e., 


&nm  — 


Cnm 


1 4-rnm(0O)#o) 


Onm  |2  (1-  |  rnm(6>o,^(})  I2) 


(2.91) 


To  estimate  mismatch  losses  in  spherical  arrays  and  to  compare  them  with  an  equivalent 
planar  array,  we  assume  a  typical  practical  case  where  the  elements  are  perfectly  matched  at 
broadside  and  have  a  reflection  coefficient  of  0.33  at  60°  off-broadside  scan.  This  corresponds 
to  a  VSWR  of  2  :  1  at  0O  =  60°.  For  the  in-between  scan  angles  we  approximate  |  Tnm  \2= 
0.128  sin  an. 

For  the  planar  array  with  Tnm  =  T  and  for  maximum  gain,  i.e., uniform  excitation  (cnm  = 
1),  the  mismatch  losses  are 

Lp  =  1-  |  Wo  =  60°)  |2  =  1125 


which  is  about  0.51  dB. 

To  evaluate  (2.91)  for  a  spherical  array,  we  replace  J2nm  by  fo"  fomax  R2  sin  a  da  df3  which 
for  uniform  excitation  gives 


jamax  |  |2  gjn  a(j[a 

Jq™*  |  a(a )  |2  (1—  [  T(a)  |2)sinada 


(2.93) 


i 


l 


W 
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where 


“W  =  TTfR'  (2'94) 

Specifically,  for  "equal  amplitude  excitation”,  c(a)  =  1,  from  (2.93)  we  get  L  =  0.345 
dB  while  if  we  apply  maximum  gain  excitation,  i.e.,  c(a)  =  cos1/2  a,  the  mismatch  losses, 
L  =  0.322  dB  are  almost  half  those  of  the  equivalent  planar  array. 


2.5  Numerical  Simulation 

2.5.1  Computer  Code 

A  Matlab  computer  code  was  generated  to  evaluate  the  radiation  characteristics  of  a  spherical 
array  and,  for  comparison,  an  equivalent  planar  array.  The  equivalent  planar  array  has  the 
same  circular  aperture  size  as  the  projected  active  sector  on  the  z=constant  plane  of  the 
spherical  array.  The  planar  array  has  a  uniform  triangular  grid  with  the  same  average  cell 
area  as  that  of  the  spherical  array.  The  code  calculates  the  radiation  patterns  in  the  0=const 
plane  vs.  elevation  angle  6  as  a  function  of  geometry,  element  type,  and  excitation.  We 
selected  a  typical  element  arrangement  on  a  sphere  consisting  of  a  number  of  rings  equally 
spaced  in  6.  Array  elements  are  then  uniformly  distributed  along  each  ring  so  that  the 
element  spacing  in  <f>  is  as  close  as  it  can  be  to  half  of  the  free  space  wavelength  (A0/2).  Each 
ring  of  course  contains  a  different  number  of  elements,  and  the  periodic  element  spacing 
on  each  ring  varies  slightly  from  ring-to-ring.  Each  element  is  designated  by  two  indeces 
( nm )  and  its  position  on  the  sphere  by  two  angles,  (an,Pnm)  as  illustrated  in  Figure  B.2  and 
Figure  B.3.  The  code  can  be  easily  adapted  to  any  array  element  distribution  on  the  sphere 
by  setting  new  values  for  (an,/3nm). 

The  radiation  elements  are  specified  by  the  element  pattern,  i.e.,  the  radiated  field  of  a 
single  element  in  a  mutually  coupled  spherical  array  environment.  In  the  code  the  user  can 
select  one  of  three  typical  choices:  isotropic,  or  element  patterns  as  given  by  either  (2.22) 
or  (2.55).  The  code  can  easily  adapt  to  any  other  element  pattern  values  whether  they  are 
theoretical,  computed  by  other  codes  or  measured  data. 

There  are  two  choices  for  the  excitation,  uniform  or  maximum  gain  excitation,  see  Ap¬ 
pendix  C.  Thus,  since  there  are  three  element  patterns,  we  have  three  different  maximum 
gain  excitations.  In  addition  to  these  preset  choices,  the  user  can  input  his  own  excitation  if 
a  low  sidelobe  level  is  desired  or  for  adaptive  pattern  control. 

In  addition  the  user  can  select  either  left  or  right  CP  or  EP,  the  beam  pointing  angles, 
and  the  size  of  the  active  sector. 

The  excitation  of  the  equivalent  planar  array  is  unity. 

The  Matlab  file  is  ”SphereArray6.m”.  The  code  input  parameters  are: 
f0=center  frequency  (GHz) 
f=operating  frequency  (GHz) 

R=sphere  radius  (m) 

ActAngleDeg=active  sector  angle  amax  (deg)  -  typical  60°. 

Th0=main  beam  direction  in  elevation  (deg) 

Phi0=main  beam  direction  in  azimuth  (deg) 


PhiDegCut=  (<?= const  ant)  plane  in  which  the  pattern  will  be  plotted 
dth=elevation  ring  spacing  (m) 
dphi=azimuth  element  spacing  in  a  ring  (m) 
sign= 

1,  for  LHCP  and  LHEP 

-1,  for  EHCP  and  RHEP 
Polar= 

1,  for  circular  polarization 

0,  for  elliptical  polarization 
Excite= 

1,  for  unity  excitation,  |  c„m  |=  1. 

2,  for  I  Cnm  |=  V cos  8, 

3,  for  I  Cnm  1=  yj (cos  $^cos2  5  +  l)/2 
Element= 

1,  for  isotropic  element 

2,  for  element  as  given  by  element  pattern  relation  (2.22),  gnm  =  cos  8nm 

3,  for  element  as  given  by  element  pattern  relation  (2.55),  gnm  =  cos  8nm(cos2  0nm+l)/2 
MissLoss= 

0,  no  mismatch  losses _ 

1,  mismatch  losses,  yl—  |  Ta  I2) 

The  output  of  the  code  is  a  plot  showing  the  array  gain  patterns  of  both  spherical  as  well 
as  planar  arrays. 

2.5.2  Calculated  Array  Patterns 

Below  we  use  the  code  described  in  the  previous  Subsection  to  demonstrate  some  of  the 
features  of  spherical  arrays.  In  all  illustrations,  sphere  radius,  R  =  10m.  The  array  grid, 
shown  in  Figure  B.2  and  Figure  B.3,  as  described  above  consists  of  number  of  array  rings. 
The  z=constant  rings  are  periodic  in  6  with  spatial  period  RdQ  =  Ao/2.  The  elements  on 
each  ring  are  also  uniformly  distributed  in  (p  however,  the  spatial  period  varies  slightly  from 
ring  to  ring  with  an  average  period  of  Esin  and<p  ~  Ao/2.  The  number  of  elements  on  each 
ring  is  determined  so  that  the  element  spacing  in  the  direction  is  as  close  to  A0/2  as  possible 
while  maintaining  2tt  periodicity.  Thus,  the  average  cell  area  is  (Ao/2  x  Ao/2)  resulting  in 
the  same  element  area  gain  Gcen  —  tt. 

The  center  frequency  is  /o  =  2  GHz  in  all  cases  unless  otherwise  stated.  The  active  sector 
is  defined  by  the  cone  angle  amax  =  60°. 

Figure  2.44  to  Figure  2.47  compare  the  radiation  pattern  of  the  above  described  spherical 
array  with  an  equivalent  planar  array  that  has  the  same  circular  aperture  as  that  of  the  pro¬ 
jected  active  spherical  sector.  Thus,  the  radius  of  the  equivalent  planar  array  is  i?sinam!22. 
Both  arrays  have  identical  elements  and  the  same  cell  area,  i.e.  the  same  element  gain 
Gcdl  =  7T- 

Specifically,  Figure  2.44  compares  the  radiation  patterns  of  a  spherical  array  with  an 
equivalent  planar  array  where  the  element  pattern  is  CP  as  given  by  (2.22).  Consequently, 
the  element  gain  for  both  arrays  is  gnm  =  Gceu  cos  8nm  =  ttcos  8nm.  Furthermore  both  arrays 
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Figure  2.1:  Schematic  view  of  spherical  array 
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Figure  2.2:  Planar  array  with  circular  aperture 
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Figure  2.4:  Spherical  array  area  normalized  to  an  equivalent  planar  array  vs.  scan  angle  9q 
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Figure  2.6:  Gain  pattern  of  96-element  circular  array 
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Figure  2.11:  Array  gain  pattern  of  triangular  array,  single  face  active, 
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Figure  2.20:  Gain  pattern  of  pentagon  array 
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Figure  2.26:  Gain  pattern  for  of  12-gon  array- 
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Figure  2.32:  Gain  pattern  of  48-gon  array 
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Number  of  sectors 


Figure  2.34:  Minimum  (blue)  and  maximum  (red)  array  peak  gain  vs.  number  of  polygon 
faces 
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Figure  2.36:  Spherical  array  showing  active  sector  in  red 
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Figure  2.37:  Active  array  area  vs.  amax 
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Figure  2,39:  Polarization  ellipse,  for  (0, 0)  element  (red)  and  (38, 0)  element  (blue) 


Figure  2.40:  Polarization  ellipse,  for  (0, 0)  element  (red)  and  (76, 0)  element  (blue) 
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Figure  2.41:  Polarization  ellipse,  for  (0,0)  element  (red)  and  (38,16)  element  (blue) 
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Figure  2.42:  Polarization  ellipse,  for  (0,0)  element  (red)  and  (76,29)  element  (blue) 
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Radiation  Pattern 


Figure  2.47:  Radiation  Pattern,  0  =  0°  cut,  Ci  =  Vcos ?,  expanded  scale 


Radiation  Pattern 


Figure  2.48:  Radiation  Pattern,  9q,  tj>o)  =  (60°,  45°),  <p  —  45°  cut 
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Radiation  Pattern 


Angle,  9  (deg) 

Figure  2.49:  Radiation  Pattern,  </>  =  0°  cut 
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Figure  2.50:  Radiation  Pattern,  <f>  =  60°  cut 
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Chapter  3 


Geodesic  Dome  Phased  Array 
Antenna 

3.1  Antenna  Architecture 

In  the  previous  chapter  we  showed  that  a  spherical  array  is  the  optimum  geometry  for 
hemispherical  coverage  in  terms  of  number  of  elements,  which  directly  translates  into  array 
cost,  and  performance.  With  all  these  advantages,  however,  spherical  array  have  not  been 
used  in  practice  primarily  because  the  fabrication  and  assembly  of  the  Beamforming  Network 
(BFN),  are  much  more  difficult  than  for  architectures  based  on  planar  array  geometry. 

In  this  chapter  we  propose  a  design  architecture,  the  Geodesic  Sphere/Dome  (GS/D) 
Phased  Array  Antenna,  that  preserves  all  the  advantages  of  spherical  arrays  while  its  fabri¬ 
cation  is  based  on  well  developed,  easily  manufactured  planar  array  technology  [2,  3].  The 
GD  phased  array  antenna  structure  consists  of  a  number  of  near-equilateral  triangular  pla¬ 
nar  subarrays  arranged  in  an  icosahedral  geodesic  dome  configuration.  It  is  important  to 
mention  that  the  Geodesic  Sphere  Phased  Array  antenna  architecture  is  not  limited  to  the 
proposed  icosahedron  geometry  ,  but  can  be  derived  from  any  one  of  the  five  regular  poly- 
hedra  (platonic  solids)  or  one  of  the  15  semi-regular  polyhedra  (Archimedian  solids).  As 
an  illustrative  example,  a  geodesic  dome  phased  array  antenna  is  shown  in  Figure  3.1.  It 
consists  of  300  near-equitriangular  flat,  panel  subarrys  obtained  by  4-frequency  subdivision 
of  the  icosahedron.  The  dome  rests  on  a  building  housing  a  computer  and  other  operational 
and  support  equipment. 

This  ’’faceted”  dome  antenna  provides  full  hemispherical  coverage  for  gateway  and  control 
(TT&C)  in  satellite  communication  systems  and/or  radar  for  air/space  surveillance  while 
at  the  same  time  exhibiting  the  following  advantages  over  the  existing  pyramid-like  and 
conformal  structures: 

1.  Keeping  the  architecture  locally  planar  and  globally  spherical  allows  the  array  to  be 
constructed  in  flat  panel  pieces  and  then  assembled  to  form  the  hemispherical  dome  structure. 
This  is  technically  and  economically  viable  because: 

-  the  design  overcomes  the  fabrication  complexity  associated  with  conformal  (curved) 
structures  since  the  subarray  fabrication  is  based  on  well-developed,  easily  manufactured 
planar  array  technology. 
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-  the  geodesic  dome  is  an  easily  realized  structure  with  well-known  mechanical  design 
and  fabrication  techniques. 

2.  The  Geodesic  Dome  antenna  preserves  all  the  advantages  of  spherical  phased  arrays 
[3]  such  as: 

-  uniform  beams  over  a  hemisphere 

-  high  gain 

-  high  instantaneous  bandwidth 

-  low  mismatch  and  polarization  losses 

-  low  input  power 

-  low  life  cycle  cost 

-  requiring  about  20%  fewer  radiating  elements  than  other  array  configurations. 

Thus  the  proposed  Geodesic  Dome  architecture  presents  a  highly  desirable,  but  up  to 
this  day,  impractical  spherical  array  design  that  provides  hemispherical  coverage  and  conse¬ 
quently  can  be  used  for  communications  and  control  with  multiple  satellites  over  the  entire 
sky.  This  optimum  design  is  achieved  by  marrying  two  mature  and  well  established  tech¬ 
nologies:  flat  panel  phased  array  technology  and  geodesic  dome  structures. 

Although  in  this  report  we  consider  the  (spherical)  geodesic  dome  antenna,  another 
similar  array  configuration,  the  geodesic  dome/cylinder  that  provides  hemispherical  coverage 
is  shown  in  Figure  3.2.  The  top  section  (hemisphere  or  cap)  is  the  same  as  in  Fig.  1,  while 
the  bottom  part  of  Figure  3.1  is  replaced  by  the  faceted  cylinder.  This  configuration  provides 
about  0.5  dB  more  gain  at  the  horizon  than  full  GD  configuration  which  is  highly  desirable 
in  satellite  communications.  In  addition  the  dome/cylinder  architecture  may  be  easier  to 
construct.  However,  in  contrast  to  full  dome,  off-zenith  beams  are  non-uniform. 
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3.2  Geodesic  Dome  Structure  -  Background 

As  previously  mentioned,  the  geodesic  sphere  phased  array  antenna  structure  is  based  on 
one  of  the  five  regular  polyhedra  or  on  one  of  the  15  semi-regular  polyhedra,  all  of  which 
could  be  inscribed  in  a  sphere  [8,  9,  10,  11].  The  five  regular  polyhedra,  also  known  as 
platonic  solids,  have  edges  of  the  same  length  and  identical  polygon  faces.  They  are  the 
tetrahedron,  cube,  octahedron,  dodecahedron,  and  icosahedron.  Most  geodesic  domes  are 
based  on  the  spherical  icosahedron,  which  is  the  result  of  the  stereographic  projection  of 
the  icosahedron  onto  its  circumscribed  sphere.  The  icosahedron  shown  in  Figure  3.3  has  20 
equilateral  triangular  faces,  12  apexes,  and  30  edges.  Thus,  the  use  of  the  icosahedron  as  the 
basis  for  the  geodesic  subdivision  of  a  sphere  results  in  the  division  of  the  spherical  surface 
into  20  equilateral  spherical  triangles.  This  is  the  maximum  number  of  equilateral  triangles 
into  which  a  sphere  can  be  divided.  For  a  geodesic  dome  antenna  requiring  hemispherical 
coverage,  the  icosahedron  can  be  truncated  as  shown  in  Figure  3.4.  For  most  practical 
applications,  however,  an  icosahedron  is  a  poor  or  too  crude  approximation  of  a  spherical 
surface.  The  spherical  surface  can  be  made  smoother  by  increasing  the  number  of  flat  panels 
(facets).  This  can  be  achieved  by  subdividing  each  of  the  20  equilateral  triangular  faces 
of  the  icosahedron  (also  called  icosa)  into  a  number  of  smaller  triangles,  or  subarrays  in 
antenna  terms.  The  number  of  triangles  depends  primarily  on  the  required  geodesic  dome 
array  performance  and  subarray  architecture,  and  fabrication  complexity  and  cost.  This 
subdivision  is  referred  to  as  a  frequency,  which  is  the  number  of  segments  into  which  the 
icosahedronal  edge  has  been  divided  by  a  particular  type  of  breakdown. 

There  are  two  basic  types  of  subdivision,  Class  I  (Alternate  breakdown)  and  Class  II 
(Triacon  breakdown).  Class  I  subdivision  can  be  realized  by  two  principal  methods:  Method 
1  which  is  based  on  similarity  of  triangle  shape  and  Method  2  based  on  similarity  of  triangle 
size. 

In  Class  I  Method  1,  and  a  breakdown  scheme  with  frequency  v,  each  edge  of  the  spherical 
icosahedron  is  divided  into  v  equal  segments,  and  the  points  of  subdivision  are  joined  with 
a  three-way  grid.  Each  face  of  the  spherical  icosahedron  is  thus  divided  into 

N*  =  £(2b  - 1)  ’  (3.1) 

n=l 

smaller  planar  equilateral  triangles,  Figure  3.5.  The  division  points  are  then  radially  pro¬ 
jected  (pushed)  outward  onto  the  surface  of  its  circumscribing  sphere,  see  Figure  3.6.  Con¬ 
necting  these  points  with  straight  lines  produces  20 N&  planar  triangles  (facets/subarrays) 
approximating  the  spherical  surface.  It  must  be  pointed  out  that  the  triangles  resulting  from 
such  subdivision  are  no  longer  equilateral  and  that  the  number  of  different  triangles  increases 
with  frequency  of  subdivision.  Class  I,  Method  2  is  the  same  as  Class  I,  Method  1,  except 
instead  each  edge  of  equilateral  triangle  i.e.,  icosa  is  divided  into  u  equal  segments,  here  the 
subtended  angle  by  the  edge  of  each  icosa  face  is  divided  into  u  equal  angles. 

We  illustrate  both  methods  on  a  simple  example,  icosa  with  Zv  subdivision  frequency. 
The  basic  principle  of  subdivision  for  Method  1  is  shown  in  Figure  3.7.  In  this  figure,  AB 
is  the  edge  or  side  length  of  the  equilateral  triangle  of  Figure  3.5.  The  edge  is  divided 
into  three  equal  sectors  where  AC  =  CD  =  DB  =  AB/ 3.  The  spherical  coordinates  of 
respective  vertexes  and  chord  factors  are  given  in  Table  3.1.  Similarly,  Method  2,  divides 
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the  9  coordinates  equally  down  the  edge  of  the  icosa  triangle,  making  all  edge  struts  equal, 
as  shown  in  Figure  3.8.  In  this  figure  a  =  (3  =  7  =  63.435/3.  Table  3.2  shows  the  respective 
spherical  coordinates  of  the  vertices  and  chord  lengths. 

To  observe  the  basic  characteristics  of  both  methods,  we  unfold  a  spherical  sector  corre¬ 
sponding  to  a  single  icosa.  Twenty  of  either  can  be  connected  together  to  make  a  spherical 
dome.  By  simply  observing  the  unfolded  diagrams,  Figure  3.9  and  Figure  3.10,  we  see  that 
in  the  case  of  equal  chord  factors  along  the  icosa  edge  (Method  2),  the  small  triangles  vary 
less  in  size,  but  vary  more  in  shape.  Disregarding  for  the  moment  the  corner  triangles  A, 
which  are  always  a  special  case,  we  notice  that  the  rest  of  the  icosa  face  is  shared  by  three 
each  of  two  quite  noticeably  isosceles  triangles,  B  and  C.  But  when  we  vary  the  edge  chord 
factors  (Method  1)  the  corner  triangles  D  stay  virtually  the  same  shape  as  before,  while  the 
rest  of  the  icosa  face  divides  itself  into  six  identical  triangles  E,  which  come  much  closer  to 
being  equilateral  than  either  B  or  C  in  the  other  diagram  (see  Figure  3.9). 

In  short,  by  letting  the  triangle  sizes  vary  more  we  can  keep  their  shapes  more  nearly 
constant.  This  seems  preferable  structurally  as  well  from  an  array  element  arrangement 
perspective.  To  elaborate  on  the  later,  it  is  known  that  for  hemispherical  coverage,  the 
array  grid  must  be  equilateral  triangular.  This  means  that  the  subaxrays  are  also  equilateral 
triangular.  On  the  other  hand,  for  efficient  use  of  the  aperture  area,  it  is  desirable  to  place 
the  largest  possible  subarray  into  a  subicosa  triangle  (dome  facet).  This  can  be  achieved 
when  a  subicosa  is  as  nearly  equilateral  as  possible.  In  this  case,  the  nonuniform  space 
around  the  subarray  is  minimum. 

Odd  and  even  frequencies  can  be  used  for  Class  I.  Figure  3.11  illustrates  the  subdivision 
of  a  spherical  equilateral  triangle  with  (a)  2-frequency,  (b)  4-frequency  and  (c)  6-frequency 
subdivision.  The  characteristic  of  this  subdivision  type  is  that  the  dividing  lines  are  parallel 
to  the  edges  of  the  original  triangles.  A  full  dome  using  Class  I,  Method  I,  with  subdivision 
frequency  u  =  2  is  shown  in  Figure  3.12. 

The  Triacon  breakdown  scheme  (Class  II)  is  similar  to  the  Alternate  breakdown  scheme 
except  the  dividing  fines  are  perpendicular  to  the  edges  of  the  polyhedron  face,  Figure  3.13. 
As  noted  in  [10],  the  Class  II  scheme  does  not  really  break  down  the  icosa  face  from  which 
it  derives  its  frequency  count,  but  smaller  and  nonequilateral  triangle  which  consists  of  two 
icosa  symmetry  triangles  ADE  (in  red)  and  CDE  back  to  back,  forming  a  diamond  ADEC, 
as  also  shown  in  Figure  3.13  for  2v  subdivision  frequency.  A  full  sphere  with  2v  Class  II 
subdivision  is  shown  in  Figure  3.14.  The  two  icosa  facets  are  indicated  with  broken-fines, 
and  the  Class  II  triangle  is  filled  in  red  color. 

Further  subdivision  of  the  icosa  face  is  possible  according  to  the  same  scheme  described 
above.  An  example  of  4i/  subdivision  is  shown  in  Figure  3.15. 

Since  Class  II  breakdown  is  based  on  symmetrical  subdivision,  only  even  frequencies  are 
possible. 

For  geodesic  dome  phased  array  applications  the  Alternate  breakdown  is  superior  because 
the  edges  of  the  triangles  fie  on  an  equatorial  circle,  which  leads  to  a  simple  design  of 
a  hemisphere  with  planar  connections.  The  connection  is  an  important  criterion  in  the 
design  of  a  geodesic  dome  antenna.  Planar  connection  cannot  be  obtained  with  the  Triacon 
breakdown  scheme.  Moreover,  in  the  Alternate  breakdown  the  resulting  triangles  are  more 
nearly  equilateral  than  in  Triacon  breakdown,  which  is  preferable  in  subarray  design. 
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Why  icosahedron? 

As  already  mentioned  above,  the  geodesic  dome  antenna  could,  in  principle,  be  con¬ 
structed  from  any  regular  or  semi-regular  polyhedra.  We  selected  the  icosahedron  as  the 
basis  for  the  geodesic  dome  antenna  structure  because  of  its  distinct  structural  and  elec¬ 
trical  characteristics.  In  comparison  to  other  polyhedra,  icosahedron-based  domes  have  the 
following  advantages  [10]: 

(1)  lower  subdivision  frequency  needed  for  comparable  average  chord  factor,  which  results 
in  a  smaller  number  of  required  subarrays; 

(2)  angles  uniformly  averaged  among  nearly  equilateral  triangles  and  less  variation  in  the 
lengths  of  members  for  a  comparable  average  chord  factor,  which  results  in  a  more  uniform 
size  and  shape  of  the  triangular  subarrays; 

(3)  curvature  is  more  uniformly  distributed,  preserving  uniform  beams  over  a  hemisphere; 

(4)  simple  truncation  to  a  planar  surface,  allowing  planar  support  structure. 

Which  Class? 

Class  I  keeps  component  variations  within  narrower  bounds,  but  at  the  cost  of  employing 
a  larger  inventory  of  different  strut  lengths.  It  permits  odd  as  well  as  even  frequencies.  Class 
I  icosahedra  have  uninterrupted  equators  at  all  even  frequencies  and  the  edges  of  the  triangles 
lie  on  an  equatorial  great  circle  which  facilitates  a  simple  hemisphere  connection  with  planar 
supporting  structure.  In  the  alternate  breakdown,  Class  I,  Method  1,  the  resulting  triangles 
are  more  nearly-equilateral  than  in  Class  I,  Method  2. 

Class  II  yields  even  frequencies  only.  Used  with  icosahedron,  it  introduces  some  trunca¬ 
tion  problems,  since  any  equator  breaks  through  some  of  the  subdividing  triangles. 

Which  frequency? 

The  subdivision  frequency  depends  on  acceptable  sphericity  and  the  number  and  size  of 
the  subarrays.  In  general,  to  minimize  the  component  inventory  we  use  the  lowest  frequency 
we  can.  For  our  10  m  geodesic  dome  array,  the  driving  factor  is  fabrication  of  the  largest 
possible  subarray.  Triangular  panels  of  1  m  on  a  side  seem  to  be  a  challenging  but  not 
impossible  task.  With  this  in  mind  we  found  that  subdivision  frequency  of  6  gives  very  good 
sphericity,  with  a  relatively  small  number  of  facets  (675)  of  reasonably  large  size. 

Decision! 

Based  on  the  above  factors,  for  our  10m  dome  we  selected  icosahedron  using  Class  I, 
Method  1  breakdown  scheme  with  subdivision  frequency  v=6.  The  Geodesic  Dome  phased 
array  antenna  is  shown  in  Figure  3.16.  More  practical  design,  a  hemisphere  with  cylindrical 
base  is  shown  in  Figure  3.17. 

There  are  numerous  publications  and  computer  codes  giving  all  the  details  (lengths  and 
angles)  for  different  configurations  and  subdivision  types.  A  computer  algorithm  for  obtain¬ 
ing  the  coordinates  of  vertices,  chord  factors,  and  dihedral  angles  of  ^-frequency  spherical 
icosahedral  domes  is  given  in  [10]. 
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Figure  3.2:  Geodesic  dome/cylinder,  subdivision  frequency  v 


Figure  3.5:  Equilateral  triangle  illustrating  Class  I-Method  1,  3r/  subdivision 
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Figure  3.7:  3  v  Ieosa  illustrating  Class  I-Method  1  subdivision 
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2,0 

43.358 

0 

1,0  / 1,1 
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3,0 

63.435 

0 
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0.412 

1,1 

20,077 

72 

2,1 
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36 

3,1 

59.008 

22.386 

2,2 

43.358 

72 

3,2 

59.008 

49.614 

3,3 

63.435 

72 

Table  3.1:  Coordinates  of  vertices  and  chord  factors  (R=l)  for  Icosa,  Method  1,  v= 3 
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0A=1  AB=1.05  a=P=7=63.435°/3 

A  (0,0) 


Figure  3.8:  Z  v  Icosa  illustrating  Class  I-Method  2  subdivision 
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63.435 
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Table  3.2:  Coordinates  of  vertices  and  chord  factors  (R=l)  for  Icosa,  Method  2,  u=Z 
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Figure  3.9:  3  u  Icosa,  Method  1,  Equal  chords  before  projection 


Figure  3.10:  3  v  Icosa,  Method  2,  Equal  edge  (angle)  division 
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Figure  3.12:  Geodesic  Dome,  v  =  2 
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B  C  B  C 

Figure  3.13:  Class  II-  2v  subdivision 


Figure  3.14:  Class  II-  2z/  subdivision  -  full  sphere,  from  [10 


Figure  3.15:  Class  II-  4i /  subdivision 


Figure  3.16:  Class  1-6 v  subdivision,  Sphere 
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Chapter  4 

Design  of  Geodesic  Dome  Structure 


4.1  Global  coordinate  systems 

A  global  coordinate  system  (x,y,  z)  with  unit  vectors  (x,y,z)  with  the  origin  at  the  center 
of  the  dome  0(0, 0, 0)  is  shown  in  Figure  4.1.  The  respective  spherical  coordinates  (r,  0.  <fi) 
with  unit  vectors  (r ,  0, <f>)  are  also  shown  in  the  same  figure. 

4.1.1  Counting  icosa’s 

Figure  4.2  shows  the  icosahedron  in  global  coordinates.  As  already  mentioned,  the  icosa¬ 
hedron  consists  of  20  equilateral  triangular  faces  (icosas).  It  has  20  apexes  and  30  edges. 
The  icosas  are  arranged  in  three  rows  designated  by  index  K  counted  from  top  to  bottom. 
First  row,  (K  =  1)  also  called  top  cap,  contains  5  icosas.  Index  L  specifies  the  location  of 
icosa  in  Kth.  row  in  +<p  (CCW)  direction  starting  at  x-z  plane.  The  middle  section  of  the 
icosahedron  ( K  =  2)  has  10  icosas,  while  the  bottom  cap  ( K  —  3)  also  contains  5  icosas.  In 
short,  each  icosa  is  designated  by  two  indeees  ( K ,  L)  where  for  K  =  1  and  K  =  3,  L  =  1  to 
5  and  for  K  =  2,  L  =  1  to  10,  as  shown  in  Figure  4.3  and  Figure  4.4. 

4.1.2  Counting  sub-icosa’s 

For  6i /  subdivision,  there  are  36  subicosa’s  in  each  icosa,  Figure  4.5.  The  vertices  are  des¬ 
ignated  by  a  two  digit  marker  as  shown  in  Figure  4.5,  the  top  vertex  being  (0,0).  The  last 
point  in  the  last  row  bottom  row  is  (6,6).  The  icosas  are  also  designated  by  two  indices 
(N.  M)  where  N  denotes  the  row  and  M  denotes  the  column  starting  at  the  tip  of  the  icosa, 
see  Figure  4.6.  Note  that  for  icosas  facing  up  the  (N,  M)  runs  from  top  to  bottom,  left  to 
right  whereas  for  icosas  facing  down  ( N ,  M)  runs  from  bottom  to  top  and  right  to  left.  In 
this  manner  the  ( N ,  M)- th  sub-icosas  will  be  identical  in  all  20  icosas. 

4.1.3  Class  I5  Method  15  6t^5  icosahedral  coordinates 

We  need  chord  factors  only  for  the  breakdown  edges  that  lie  partially  or  wholly  inside  the 
symmetry  triangle.  For  6 v  subdivision  this  means  only  the  twelve  edges  drawn  in  red  in 
Figure  4.5  for  A#l,  thus  twelve  out  of  sixty-three.  Thus,  the  only  coordinates  we  need  are 
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the  end  points  of  those  twelve  edges;  a  total  of  ten  coordinates.  The  rest  of  the  system  repeats 
symmetrically;  thus  (0, 0)  -  (1, 0)  corresponds  to  (5, 0)  -  (6, 0);  (6, 0)  -  (6, 1)  corresponds  to 
(6, 5)  —  (6, 6)  etc.  For  icosahedron,  Class  I,  Method  1  subdivision,  the  spherical  coordinates 
of  these  vertices  are  given  in  Table  4.1  The  rest  of  the  vertex  coordinates  in  A  #1  can 
be  deduced  from  the  symmetry  consideration  as  mentioned  above.  For  convenience  they 
are  listed  in  Table  4.2.  Coordinates  of  triangle  (icosa)  No.  2  shown  in  Figure  4.7  can  be 
generated  from  data  of  triangle  No.  1  by  rotation.  They  are  also  given  in  Table  4.2. 

With  known  coordinates  of  triangles  #1  and  #2,  triangle  sections  3,4,5  and  6  which  are 
marked  in  red  in  Figure  4.8  and  Figure  4.9  can  be  similarly  obtained  according  to 
Z:<f>  =  fa  +  72°,  6  =  0i 
4:  4>  =  fa  -  36°,  9  =  180°  -  02 

5:  (j)  —  fa  +  36°,  9  =  same  as  A#4 

6:  0  =  72°  +  fa,  0  =  02- 

From  this,  the  rest  of  the  coordinates  of  the  icosahedron  can  be  obtained  by  rotation 
using  symmetry  of  the  structure.  Thus,  each  triangle  is  specified  in  space  by  three  points 
(comer  coordinates)  Pi{R,  9i,  fa),  Pz{R,  02,  fa)  and  P$(R,  9$,  fa). 

Once  the  spherical  coordinates  ( R ,  0,  fa  of  each  triangle  are  known,  the  respective  co¬ 
ordinates  in  a  Cartesian  system  P\{x\,y\,zi),  P2(x2ty2,z2)  and  P${xz,yz,zz)  can  be  found 
from 

x  =  .Rsin0cos0 

y  —  i?sin0sin0 

z  =  Rsin9.  (4.1) 

4.1.4  Class  I,  Method  1,  6 z/,  chord  lengths 

The  chord  factors  can  be  now  evaluated  from 

d  =  \J2  —  2  [cos  9i  cos  d2  +  cos{fa  —  fa)  sin  9i  sin  $2]  (4.2) 

where  for  0,  <f>  we  insert  values  from  the  above  table.  Thus  to  calculate  the  distance  from 
(0, 0)  —  (1, 1)  we  use  0  for  values  of  fa,  fa  and  9i,  and  9.3247035  for  02.  The  chord  factor 
turns  out  to  be  0.162567.  Following  the  same  procedure,  the  rest  of  the  chord  factors  can 
be  evaluated.  They  are  given  in  Figure  4.10  for  a  sphere  radius  R  =  1.  For  any  other  R, 
the  chord  lengths  must  be  multiplied  by  R.  The  chord  factors  (R  =  1)  or  the  chord  lengths 
(R  >  1)  can  be  also  evaluated  from  following  equations: 

di:j  =  yj ( x{  -  Xj)2  +  (j/i  -  7/j)2  +  {zi  -  Zj)2  (4.3) 

where  as  already  mention  (x,,yi,Zi),  i  —  1,2,3  and  ( Xj,yj,zj ),  j  =  1,2,3  are  coordinates  of 
triangles. 

4.1.5  Sub-icosas  (panel)  types 

For  61/  subdivision,  as  illustrated  in  Figure  4.11,  there  are  6  different  panel  types,  which  for 
convenience  we  designated  them  from  type  1  to  type  6.  Figure  4.11  displays  all  triangle  types 
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with  different  color.  Notice  that  there  are  two  asymmetrical  (mirror)  triangles  for  types  3 
and  5.  Triangle  types  1,2  4  and  6  lie  on  the  icosa’s  symmetry  axis,  and  because  of  that  they 
are  isosceles  triangles.  The  smallest  triangles  are  ’’corner”  triangles  (dark  blue)  while  the 
largest  triangle  lie  in  the  middle  and  are  marked  in  red. 

4.2  Local  coordinate  system 

We  introduce  a  local  coordinate  system  (%nm,  Vnm,  z^m)  with  respective  unit  vectors 
(*NM,yNM,ZNM)  centered  at  ONm{Rnm,olnm,0nm)  or  0NM{xQ,  y0,zo),  Figure  4.12.  The 
local  coordinate  system  is  oriented  so  that  xjvm  =  6,  Ynm  =  4>,  %nm  —  r  at  Omm- 


4.2.1  Sub-icosa  (panel)  in  local  coordinate  system 

We  consider  the  (JV,  M) th  subicosa  or  triangle  (panel)  with  coordinates  Pi(xi,yi,  zy), 

P2(x2.  y2,  z2)  and  P3(x3,  y3, 23).  We  define  the  unit  normal  to  the  plane  n  =  xnx  +  yny  + 
znz]  directional  cosines  (nx,  ny,n2)  to  be  determined.  Then,  we  define  the  line  LL,  with 
the  same  directional  cosines  (nx,ny,nz)  passing  through  the  global  origin  O.  The  origin 
OpfM(xQ,y0,zQ)  of  the  local  coordinates  is  specified  by  the  intersection  of  line  LL  and  the 
triangle  plane  as  shown  in  Figure  4.13.  The  zNM-axis  of  the  local  coordinate  system  coincides 
with  the  normal  n  while  axis  x^m,  Vnm  lies  in  the  triangle  plane. 

First  we  determine  the  directional  cosines  of  the  n  then,  the  location  of  the  Omm  i.e., 
coordinates  (xq,  y0, 20)  and  finally  the  directional  cosines  of  the  (xnm,  Vnm ,  *nm)  local  axis. 


Directional  cosines  of  normal 

The  general  equation  of  a  plane  in  space  is 

Ax  +  By  +  Cz  +  D  =  0.  (4.4) 


Specifically,  for  the  plane  passing  through  three  points  in  space  Pi(xi,yi,zi),  ^2(^2, 2/2,  z2), 
P$(%3, 2/3,  z3)  we  may  write 


or 


x-xi  -  y  —  yi 

%2  ~  xi  2/2  —  2/1 

-  x 1  y3  —  yi 


2  —  2x  • 
Z2-Z1 
Z3  -  z 1 


=  0. 


(x  -  Xi)[(y2  -  yi)(z3  -  z{)  -  (y3  -  yx)(z2  -  zi)] 

+  (y  -  2/l )[(  0^2  -  ®l)(23  -  Zi)  +  (x3  -.Xi)(z2  -  2x)] 

+  (2  -  Zi)[(x2  -  Xx)(y3  -  yi)  -  (x3  -  xx)(y2  -  yi)]  =  0. 


In  (4.4),  we  see  that 

A  =  (y2  —  yi)(z3  -  zi)  -  (y3  -  yx)(22  -  zl) 

B  =  -[(-(x2-xi)(23-21)  +  (x3-x1)(z2-21)] 
C  =  (x2 -Xx)(y3 -yi)  -  (x3 -xx)(y2 -yi) 


(4.5) 


(4.6) 


(4.7) 
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and 


(4.8) 


D  =  —x\A  —  y\B  +  z\C. 

Thus  the  directional  cosines  of  the  normal  vector  n  to  the  plane  are: - 

A 

nx  -  cos a~  ^WTWTC^ 
ny  =  cos(3  =  -j======= 

C 

n,  =  cos  7  =  .  . 

’  VA 2  +  B2  +  C2 

Warning:  (a,/?,  7)  in  (4.9)  not  to  be  confused  with  {chnm,Pnm)  in  Figure  4.13. 


ny  =  cos  (3  = 


nz  =  cos  7  = 


Location  of  Onm' 

The  general  equation  of  a  line  in  parametric  form  passing  through  global  origin  O  is 

x  —  nxt 
.  y  =  nyt 
z  —  nzt. 

Since  OAw(ro,  2/0,  zQ)  lies  on  this  line  we  can  write 

Xq  = 

2/0  =  nyto 

ZQ  =  nzt0.  -- 

At  the  same  time  Onm  also  lies  on  the  planar  triangle  so  we  can  write 

Ax  0  +  Byo  +  Czq  +  D  =  0  . 

or  using  (4.11) 


From  here 


An*xto  "I-  Bnyt0  "I-  cnzt0  *4*  B  —  0. 


•  Anx  +  Bny  +  Cny ' 


Substituting  (4.14)  into  (4.11),  the  global  coordinates  of  the  origin  Onm  are 

Xq  —  nxtQ 
2/o  =  nyto 
zq  -  nzt0. 

In  spherical  coordinates  the  coordinates  of  Onm{Rq,  ocnm/Pnm)  are: 

f?~ - .  frr^  _]_  rt  ,2  I  -2 


•Ro  =  yJxl  +  yl  +  zl 


aNM  —  arccos  - 
Rq 

Pnm  =  arctan 

Rq 


(4.10) 


4i.ll>..  - 


(4.12) 


(4.13) 


(4.14) 


(4.15) 


(4.16) 

(4.17) 

(4.18) 
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Local  axis 

Consequently,  the  directional  cosines  of  local  z^Af-axis  are: 

nx(xj vm)  =  cos  a jvm  cos  (3^m 
%  {%nm)  —  cos  oiNM  sin  j3j vm 
nz(xNM)  =  —  sino/vM 

where  the  identity 


\fnl{xNM)  +  n2y{xNM)  +  n2z{xxM)  =  1 
holds.  Similarly,  directional  cosines  of  the  local  y^-axis  axe: 

nx(VNM )  =  —  sin^jvM 
ny{VNM)  =  cos  /S.vm 
^(s/wm)  =  0 

with  _ 


(4.19) 

(4.20) 


\Zni(z/.vM)  +  ny(yNM)  +  «i(yiw)  =  1 
As  mentioned  n  coincides  with  the  unit  vector  i>iM)  so 

^x(^2Vm)  =  'H'z 

Tly  {ZNM )  =  Tty 


and 


nz{zNM)  —  Uj.. 


(4.21) 

(4.22) 

(4.23) 

(4.24) 


XnM 

-M 

’  X  —  Xo  ' 

Vnm'- 

y-yb 

.  *NM  . 

N 

1 

£ 

Coordinate  transformation  -  going  from  global  to  local 

Any  given  point  in  global  coordinates  (x,y;z)  can  be  transformed  into  local  coordinates 
(xnm,Vnm,znm)  by  following  translation  -  rotation  matrix: 


(4.25) 


where 

Kx(xnm)  TIz{xhm) 

[A]  =  nx{yNM )  ny{yNM)  nz{yNM)  •  (4.26) 

.  nx(z?iM)  ny(zNM)  nz(zffM) 

Each  subicosa  triangle  can  now  be  represented  in  local  coordinates  by  using  (4.26)  for  each 
the  three  comer  points.  Because  the  unit  vector  is  perpendicular  to  subicosa,  z^m  =  0. 
Thus,  one  can  write  for  the  coordinates  of  the  subicosa  in  local  coordinate  system 

xnm  =  tt-x(xnm)(x  ~  xo)  +  ny{xhiM)(y  ~  Vo)  +  nz{xNM){z  —  z0) 

Vnm  =  nx(yNM)(x  —  xo) +  ny(yivM)(y —  yo) +  nz(yNM)(z  —  Zo) 

zmm  =  0.  (4.27) 
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Figure  4.14  to  Figure  4.29  show  all  eight  types  of  subicosas  in  local  coordinates.  As  an 
example  we  selected  the  subicosa  triangles  ( N ,  M)  to  be:  (1, 1),  (2, 1),  (2, 2),  (2, 3),  (3, 1), 
(3, 2),  (3, 3),  (4, 4).  The  specific  location  of  these  triangles  in  the  icosa  triangle  is  also  illus¬ 
trated  with  fill-in  red  color  in  Figure  4.14  to  Figure  4.29. 

For  each  triangle  point  (xnm,Vnm)  coordinates  are  also  given  in  Figure  4.14  to  Fig¬ 
ure  4.29.  In  addition  chord  lengths  (triangle  sides)  are  denoted  with  lower  case  Roman 
letters  according  to  Figure  4.10.  Chord  lengths  may  be  computed  from 

dij  —  V  (xpfMi  —  XMMj)2  +  (yNMi  —  VNMj)2  (4.28) 

where  (x^mu yNMi),  i  =  1,2,3  and  ( XNMj,VNMj ),  3  =  1,2,3  are  the  coordinates  of  the 
subicosa  triangles  in  local  coordinate  system.  These  chord  lengths  must  be  the  same  as 
those  computed  by  (4.3)  in  global  coordinates,  which  may  be  a  good  check  for  coordinate 
transformation  and  computer  code. 

For  convenience  we  determine  the  center  of  mass  point  in  local  coordinates: 

c  XffMl  +  XNM2  +  XNM3  fAO(\\ 

XNM  ~  —  o 


c  _  VNM1  +  UNM2  +  VNM3 

Vnm  ~  5 


(4.30) 


Center  of  mass  point  of  every  subicosa  is  also  shown  in  Figure  4.14  to  Figure  4.29  by  a  plus 
sign. 


4.2.2  Sub-icosa  base 

Because  the  subarray  must  be  equilateral  triangular,  here  we  determine  the  orientation  of  the 
largest  possible  equilateral  triangle  that  fits  into  a  subicosa  triangle.  Figure  4.30  shows  the 
triangle  with  sides  a,  b  and  c  in  local  coordinates.  As  illustrated  in  this  figure,  it  is  possible 
to  fit  three  typical  equilateral  triangles  within  a  subicosa  triangle  so  that  one  side  of  the 
equilateral  triangle  coincides  with  one  of  the  sides  of  the  subicosa  triangle.  However  these 
triangles  are  not  the  same.  The  largest  triangle  is  the  one  that  coincides  with  the  subicosa 
base  closest  to  the  isosceles  triangle;  side  b  in  example  of  Figure  4.30.  If  the  subicosa  is 
isosceles  triangular,  which  axe  all  triangles  along  the  axis  of  symmetry,  i.e.,  triangle  types 
1,2,4  and  6,  then  the  base  of  the  isosceles  icosa  is  also  a  side  of  the  equilateral  triangle. 
For  these  triangles  the  bases  are  sides  b,  6,  h,  and  h,  respectively.  For  triangle  type  3a  and 
3b,  the  base  is  chord  length  e,  while  for  triangle  types  5a  and  5b,  the  base  chord  length 
is  g.  Thus  each  subicosa  in  a  geodesic  dome  has  one  side  designated  as  the  base.  The 
bases  of  subicosas  (panels)  in  a  single  icosa  triangle  are  designated  in  red  in  Figure  4.31. 
For  convenience,  the  base  is  also  marked  in  Figure  4.30  by  the  line  OnmP  perpendicular 
to  the  base.  Furthermore,  using  the  same  rule,  the  bases  are  marked  for  each  subicosa  in 
Figure  4.11,  and  in  Figure  4.15  to  Figure  4.29. 
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4.3  Reference  coordinate  system 

We  define  one  more  coordinate  system,  the  reference  coordinate  system  {x'NM,y'NM,z'NM) 
shown  in  Figure  4.32  and  Figure  4.33.  A  subarray  is  designed  and  tested  in  the  reference 
coordinate  system.  The  origin  of  this  coordinate  system  is  at  Omm  while  the  z'NM  axis 
coincides  with  the  zmm  axis  of  the  local  coordinate  system.  As  illustrated  in  Figure  4.33  the 
x'NM  axis  coincides  with  line  UnmP  and  is  pointing  toward  the  tip  of  the  equilateral  triangle. 
Thus,  the  reference  coordinate  system  can  be  obtain  by  rotating  the  local  coordinates  about 
the  zmm  axis  by  angle  $mm  hi  the  CCW  direction.  For  subicosas  along  the  vertical  axis 
of  symmetry,  i.e.,  (1,1),  (2,2),  (3,3),  (4,4),  (5,5)  and  (6,6)  the  rotation  angle  is  0  or  180 
deg.  For  other  subicosas  the  rotation  angle  depends  not  only  on  azimuthal  position  but 
also  on  elevation.  Thus,  the  rotation  angles  depend  on  the  subicosas  position  within  the 
icosa  and  on  the  icosas  elevation  position.  We  compute  the  rotation  angles  with  the  help 
of  Figure  4.33,  An  equilateral  triangle  with  side  as  is  shown.  For  the  sake  of  analysis  we 
assume  that  {xNm,  yNm){xNM3,  Vnmz)  is  the  base.  The  angle  that  the  line  parallel  with 
the  base  makes  with  respect  to  the  -tx  mm  axis  is 

4%}!*  =  arctan  (4.31) 

XNMZ  -  XNM2 

Then  the  rotation  angle 

<?nm=$T  +  ’ r/2.  (4.32) 

Figure  4.35  shows  the  four  basic  icosas  and  respective  subicosas.  Each  subicosa  in  these 
four  basic  icosas  has  a  different  angle  L,  N,  M )  of  rotation.  These  angles  for  the 

four  basic  icosas  are  given  in  Table  4.3  to  Table  4.6.  These  four  icosas  constitute  the  basic 
cell  of  the  full  dome,  i.e.,  by  rotating  Figure  4.35  in  steps  of  72  deg  about  the  2-axis,  we 
would  get  full  dome.  Thus,  the  angles  <pmm  of  the  rest  of  the  subicosas  in  the  dome  can 
be  determined  from  those  in  Table  4.3  to  Table  4.6,  by  taking  into  account  periodicity  in 
azimuth.  Thus,  we  use  the  following  identities: 


<Pnm(  1,  L,N,M)  =  <pnm0-,  1)  N,  M), 

II 

to 

..,5 

(4.33) 

V>nm(%  L,  N ,  M)  =  <pmm( 2j  1,  N ,  M), 

L  =  1,3,. 

:.,9 

(4.34) 

<Pnm( 2,  £»  N,  M)  =  ipN M(2, 2,  N,  M),  *• 

1  =  2,4,. 

..,10 

(4.35) 

<Pnm(  3,  L ,  N,  M)  =  <pnm(  3, 1  ,N,M), 

t-» 

II 

to 

...5. 

(4.36) 

Thus  the  reference  coordinates  (x'NM,ytNM)  can  be  obtain  by  rotating  the  local  coordi¬ 
nates  (xnm,Vnm)  about  zNM  axis  by  angle  pNM  in  CCW  direction.  In  reference  to  Fig¬ 
ure  4.34,  this  transformation  can  be  described  by  the  following  rotation  matrix  relation: 

x'nm  ^  _  (  cos  pnm  sin  pmm  \  /  xmm 

Vnm  )  \  -  sin  pnm  cos  ipmm  )  \  Vnm 

Taking  the  inverse  of  the  rotation  matrix  the  reference  coordinates  can  be  transfered  back 
to  local  coordinates  by  using 

xnm  \  _  /  cos p mm  —sin pmm  \  (  x'nm 

Vnm  )  \  sin  pm m  cos  pnm  )  \  Vnm 


(4.37) 
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4.3.1  Sub-icosas  in  reference  coordinates 


Now,  we  transform  all  subicosa  triangles  into  reference  coordinates  according  to  the  above 
prescription,  which  for  the  three  triangle  points  becomes _ _ _ _ _ _ _ 


x'nmi  ~  xnmi  cos  <Pnm  +  y'N M1  sin  <j>nm 

VnMI  =  ~XNM1  SintfNM  +  Dn Ml  C0S  VNM  (4.39) 

XNM2  =  XNM2  COS  tyNM  +  VnM2  <PNM 

v'nm2  ~  xnm2  sin  (pNM  +  y'NM2 cos  Pnm  (4.40) 


x’nm3  =  xnm3  cos  <pNM  +  y'Nm  sin  ippjM 

Dnmz  =  — xnmz  sin  <Pnm  +  Vnmz  cos  ^Pnm-  (4-41) 

As  an  illustrative  example,  Figure  4.38  to  Figure  4.45  show  all  six  triangle  types  belonging 
to  icosa  ( K ,  L)  =  (1, 1)  in  reference  coordinates.  Numerical  values  for  the  triangle  coordinates 
are  also  given  in  these  figures.  We  see  that  all  six  types  differ  considerably  in  size.  For 
convenience  the  six  basic  types  are  plotted  on  the  same  axis  (x'NMl  y'NM)  in  Figure  4.36. 
Moreover,  for  easy  comparison,  we  plotted  the  same  triangles  so  that  the  bases  coincide  on 
the  same  horizontal  fine,  Figure  4.37.  To  efficiently  populate  these  triangles  with  radiating 
elements,  it  is  obvious  that  each  of  them  will  require  a  different  number  of  elements.  More 
on  this  subject  will  be  said  in  the  next  section. 


4.3.2  Fitting  equilateral  triangle  into  sub-icosa  triangle 

In  this  section  we  determine  the  maximum  side  length  of  the  equilateral  triangle  that  fits 
into  each  subicosa  type  and  the  respective  coordinates  of  the  equilateral  triangles  in  the 
reference  coordinate  system. 

Type  1: 

Figure  4.38  shows  the  subicosa  in  reference  coordinates.  We  see  that  the  maximum  side 
of  the  equilateral  triangle  ag  is  limited  by  the  height  of  the  subicosa.  We  find  ag  as  follows: 

b 

cos  a  =  —  (4.42) 


Frome  here 


V3  . 

ti a  —  cle~t~  =  asm  a. 


2  b 

cle  —  -~7=asin(arccos  — )  —  0.1521. 
V  3  2a 


(4.43) 

(4.44) 
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Consequently,  the  coordinates  of  the  equilateral  triangle  are: 

(XNM1E1  V'nMIe)  —  (xiVMl>0) 
(XNM2E>  VnM2e)  =  (xNM2i  qe/^) 
{^NMZE^VnMZe)  ~  ( xNM3’~aE /2)- 


(4.45) 


The  center  of  mass  is  therefore 


lEmass 


_  XNMIE  +  XNM2E  +  xNM3E 


VEmass  0* 


(4.46) 


Type  2: 

Following  the  same  procedure  using  Figure  4.39,  we  get 


and 


it  t  i) 

ce  =  -7=csin(arccos  — )  =  0.1790 
v3  2c 


(xNM1EjVnM1e)  —  (xWl)0) 

(x>NM2Ei  VnM2e)  =  ixNM2i  °f/2) 
{xNMZEjVnMZe)  =  (XNMZ1  ~aE/^)- 


(4.47) 


(4.48) 


Center  of  mass 


■'Emass 


X>NM\E  +  XNM2E  +  XNM3E 


VEmass  0. 


(4.49) 


Types  3a  and  3b:  Subicosas  of  type  3  and  5  are  the  triangles  with  all  three  sides 
different,  as  shown  in  Figure  4.40  and  Figure  4.41.  In  type  3  subicosas  the  base  is  chord 
length  e.  The  maximum  size  of  an  equilateral  triangle  is  determined  by  the  subicosas  height 
h&.  To  determine  as  we  first  find  height  h&  as  follows: 


a  =  2  arctan , 


^  (e  +  d)2  —  c2 


and  therefore 
Finally, 


h&  =  dsina. 


2  h 

aE  =  ^=  =  0.1781 


(4.50) 

(4.51) 

(4.52) 


ixNMlEi  v'nM\e) 
(xNM2E>  V'nM2e) 

{x'nmzei  Vnmze ) 


(XNMV  VnMi) 

(x'nm2->  v'nmi  +  o-e/  2) 
(p'nmziV'nmi  ~  °s/2)- 


(4.53) 
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Center  of  mass 


_  x'nm\e  +  x'nm2E  +  x 'nmze 
3 

__  VjunE  +  v'nmie  +  Vnmze 
3 


(4.54) 


Subicosa  of  type  3b  is  the  mirror  image  of  type  3a.  In  reference  to  Figure  4.41,  the 
same  expressions  as  for  a  subicosa  of  type  3a  can  be  used,  giving  the  same  a#  as  in  (4.52). 
Expressions  for  equilateral  triangle  coordinates  (4.53)  also  apply. 

Type  4: 

In  this  case  from  Figure  4.42  we  can  write 


2  h 

aE  =  -7=/sin(arccos— )  =  0.2027 
v3  2/ 


(4.55) 


Center  of  mass 


(p'nMXEiV'nMIe)  =  (XATM1)  VnMi) 

(XNM2E>  VnM2e)  =  (XNM2>  afi/2) 

{xNMZEiV'nMZe)  =  (x'nmz>  qe/2)  . 


_  XNM1E  +  x  NM2E 

~  3 

=  0. 


(4.56) 


(4.57) 


Types  5a  and  5b:  The  base  of  this  type  of  subicosas  is  chord  length  g.  In  this  case 
the  maximum  length  of  the  equilateral  triangle  is  limited  by  the  base  length,  as  seen  in 
Figure  4.43  and  Figure  4.44.  Thus,  for  triangles  of  type  5a  and  5b,  as  =  g  =  0.1980  and 


Center  of  mass 


(X>NM1E>  v'nmie)  =  (XNM11  — ' ®s/2) 

{x’nM2Ei  Vnmze)  =  (x’nmi  +  aE sin(7r/3),  0) 

(XNMZE>  VnMZe)  =  (xNMZi  aJs/2)- 


NM1E  XjjmE  +  x'nMZE 

3 

0. 


(4.58) 


(4.59) 


Isosceles  Subicosa  of  Type  6,  (see  Figure  4.45)  is  limited  by  its  base  length  h.  Thus,  in 
this  case 

aE  =  h  =  0.2154  (4.60) 


(xNMlEi  v'nmie) 
{x'nM2EiV'nM2e) 
{x'nmzei  v'nmze) 


(x'nmii  as/2) 

(xNM2i  ~aE/ 2) 

(xNMi  +  aE  sin(7r/3),0). 


(4.61) 
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Center  of  mass 


/  _  ^NMIB  +  x'nM2B  +  X>NMZE 

' amass  g 

*L»  =  0-  (4.62) 

All  six  subicosa  types  with  their  respective  equilateral  triangles  in  the  reference  coordinate 
systems  are  shown  in  Figure  4.46  to  Figure  4.61.  The  coordinates  of  the  triangles  are  also 
written  in  blue  color  for  subicosas  and  in  red  for  equilateral  triangles.  Figure  4.62  shows 
the  six  different  equilateral  triangles  in  the  reference  coordinate  system.  For  comparison 
they  are  also  shown  in  Figure  4.63  where  the  bases  are  aligned  on  the  same  line.  As  already 
mentioned  for  efficient  use  of  the  area  we  would  like  to  populate  these  triangles  with  as  many 
radiating  elements  as  possible  we  can.  On  the  other  hand,  a  different  number  of  elements 
on  each  panel  would  require  different  beamforming  networks,  which  is  undesirable  from  ease 
of  fabrication  and  low  cost  perspective.  From  Figure  4.63  it  is  seen  that  equilateral  triangle 
types  2  and  3  are  very  similar.  Also  equilateral  triangles  4  and  5  are  very  close  in  size.  As 
will  be  shown  latter,  this  means  that  we  can  reduce  total  number  of  different  triangles  from 
6  to  4. 


4.3.3  Counting  the  array  elements  in  reference  coordinate  system 

{x’nm-  v'nm) 

Similarly  as  we  counted  subicosas  in  a  single  icosa  with  two  indeces  (N,  M)  where  N  is 
counted  from  left-to-right,  M  from  top-to-bottom,  here  we  designate  each  subarray  element 
with  two  indeces  (n,  m)  where  index  m  runs  from  left-to-right,  n  runs  from  top-to-bottom,  as 
indicated  in  Figure  4.65.  The  elements  are  arranged  in  an  equilateral  grid  with  dimensions 
(<4,  dy).  There  are  nmMax  rows  in  a  subarray  (subicosa).  As  already  mentioned  we  would 
like  to  place  as  many  elements  as  we  can  in  the  equilateral  triangle  with  spacing  around  the 
edges  equal  to  equilateral  cell  size.  Considering  this, 


a>E 


nmMax  =  lower  integer(- - 1). 

dy 


(4.63) 


For  geodesic  dome  with  R  =  5m,  Table  4.7  shows  the  number  of  rows  for  each  panel 
type.  Figure  4,66  to  Figure  4.71  show  the  panel  with  radiating  elements  in  the  reference 
coordinate  system.  The  radiating  elements  are  circular  microstrip  patches. 


4.3.4  Location  of  the  array  elements  in  reference  coordinates 


Once  the  number  of  rows  is  determined  for  each  panel  (subicosa)  type  we  can  then  determine 
the  location  of  the  array  elements  in  reference  coordinate  system  {x'NM)  The  equilat¬ 

eral  array  with  nmMax  rows  is  centered  with  respect  to  center  of  mass  point  (x'massE,  y'massE) ' 
of  the  equilateral  triangle  as  shown  in  Figure  4.64.  The  radial  vector  from  the  origin  of  the 
local  coordinate  system  (0,0)  to  the  center  of  mass  point  (x'massEl  y'massE)  is 


^ NM  ^NM^massE  d~  Y NtvtV mass E  ‘ 


(4.64) 
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The  vector  from  the  center  of  mass  point  to  the  first  array  element  (ra,  m)  =  (1, 1)  is 

=  -§(x!vMtan(7r/6)  +  y'NM).  . ____(4.65) 

where  L  is  the  element  center-to-center  distance  between  the  first  and  last  elements  in  the 
first  row,  i.e.,  L  =  nmMax  dy.  Thus  the  vector  from  the  first  element  (1, 1)  to  element 
(n,m)  is 

rnm  =  +  y 'nmK™  ~  !)  +  ~  2  "R  (466) 

where  dy  is  the  equilateral  lattice  grid  side  length  and  dx  =  dy  sin(7r/3).  Therefore  the 
location  of  element  (n,  m)  in  the  reference  coordinates  is 

T'nm  —  V°NM  +  TNM  +  Tnm  (4-67) 


which  may  be  rewritten  in  the  following  form: 

v'nm)  =  X NMX'nM  +  $ 'nmVnM 

where 


and 


XNM  =  x'massE  ~  %  ^Tr/6  +  (n  -  1)4 

/  /  f  (•/  -\  (71  1)1  , 

Vnm  =  2/^K  -  2  +  l(m  -  !)  +  — 2 — 


(4.68) 

(4.69) 

(4.70) 


4.3.5  Location  of  the  array  elements  in  local  coordinates 


Element  locations  originally  given  in  the  reference  coordinates  can  be  transformed  into  local 
coordinates  according  to  inverse  rotation  matrix  relation  (4.38)  . .  - . . - . 


XNM 

COS  tfiNM 

—  sin  <pnm 

5: 

_  Vnm  _ 

_  sin  <p  jvm 

cos  ip  mm 

L  Vnm  j 

Thus,  in  local  coordinates  the  array  elements  are  located  at 


(4.71) 


*'nm(xNM,  DNm)  =  XNMXNM  +  tiNMljNM  (4.72) 

where  according  to  (4.71) 

xnm  =  x'hm  cos  Vnm  —  v'hm  sin  Vnm 

Vnm  =  x'nm  s^nVNM  +  y'MMC0SVNM-  (4-73) 
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4.3.6  Location  of  the  array  elements  in  global  coordinates 

Here  we  transfer  local  coordinates  of  the  elements  into  the  global  coordinate  system  (x,  y,  z). 
With  the  help  of  (4.25)  and  (4.26)  we  can  write""  . ...  JH.  _ _ ... 


%nm 

Vum 

-  [A]”1 

XNM 

Vnm 

4- 

Xq 

yo 

Z-nm 

.  zNM  . 

.  Z0  . 

where  [A\~l  is  inverse  of  (4.26).  Elements  of  [A]-1  are 

«X1  =  ny(yNM)nz(zNM)  ~  nz(yNM}ny(zNM) 
fl21  ~  nz{UNM  (ZNM )  iR>x{3}Nm)'R>z{zNm') 
fl3i  =  nx(yNM)ny(zNM)  -  ny(yNM)nx(zNM) 


(4.74) 


(4.75) 


fll2  —  flz  (%NM  (%NM )  Tly{%NM  )flz 

022  =  nx{xNM)nz{zNM)  —  nz(xNM)nx(zHM) 

fl32  =  Tiy{xNM)nx{zNM)  —  nx{x^M)ny{z^M)  (4-76) 


fli3  =  ny(xNM)nz(yNM)  —  nz{xNM)ny{yNM) 

023  =  nz(xNM)nx(yNM)  ~  nxi^NM^ziVNM) 

fl33  =  nx{xi<iM)ny{yNM)  -  fiy(xtfM)nx(yNM)- “ .(4  77).. 

The  radial  vector  from  the  origin  of  the  global  coordinate  system  to  the  (nm)th  element 
is  therefore 

E?im(®  i  V •  z)  ~  4”  yynm  4"  ZZnjrl .  ~ . ~  “(4.78) 
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Figure  4.1:  Global  coordinate  system  (a:,  y,  z ) 


Figure  4.2:  Icosahedron  in  global  coordinates 


Figure  4.5:  Icosa  A#1  showing  numbering  of  the  sub-icosa’s  (facets)  and  vertices 


Panel:  (N,M) 


Figure  4.6:  Icosa  showing  scheme  for  counting  the  sub-icosa’s 
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Figure  4.7:  Ieosa  No.  2  pertaining  to  evaluation  of  respective  vertex  coordinates,  Class  I, 
Method  1,  6u  frequency  [10] 


Triangle  #1,  6v 

Vertex 

0  (deg) 

$  (deg) 

0,0 

0.0000000 

0.0000000 

1,0 

9.3247035 

0.0000000 

1,1 

9.3247035 

72.0000000 

2,0 

20.0767513 

0.000Q000 

2,2 

16.4722107 

36.0000000 

2,2 

20.0767513 

72.0000000 

3,0 

31.7174744 

0.0000000 

3,1 

27.2237351 

22.3861776 

3,2 

27.2237351 

49.6138225 

3,3 

31.7174744 

72.0000000 

4,0 

43.3581976 

0.0000000 

4,1 

39.1034177 

16.0353713 

Si 

37.3773682 

36.0000000 

1 

39.1034177 

55.9646288 

KiU 

43.3581976 

72.0000000 

5,0 

54.1102453 

0.0000000 

5,1 

50.6513527 

12.4463843 

5,2 

48.4869490 

27.7323015 

5,3 

48.4869490 

44.2676986 

5,4 

50.6513527 

59.5536157 

5,5 

54.1102453 

72.0000000 

6,0 

63.4349488 

0.0000000 

6,1 

60.9162275 

10.1562304 

6,2 

59.0080312 

22.3861776 

6,3 

58.2825256 

36.0000000 

6,4 

59.0080312 

49.6138225 

6,5 

60.9162275 

61.8437696 

6,6 

63.4349488 

72.0000000 

Triangle  #2,  6v 


Vertex 


0 (deg) 

<t>  (deg) 

71.0914502 

5.7778072 

69.4828700 

16.7344818 

68.5073635 

29.3546281 

68.5073635 

42.6453719 

69.4828700 

55.2655182 

71.0914502 

66.2221928 

80.1168545 

11.8185857 

79.4443283 

23.3461239 

79.1876831 

36.0000000 

79.4443283 

48.6538761 

80.1168545 

60.1814143 

90.0000000 

18.0000000 

90.0000000 

29.8185857 

90.0000000 

42.1814143 

90.0000000 

54.0000000 

99.8831455 

24.1814143 

100.0928406 

36.0000000 

99.8831455 

47.8185858 

108.9085499 

30.2221928 

108.9085499 

41.7778072 

116.5650512 

36.0000000 

Table  4.2:  Coordinates  of  the  vertices  for  icosas  A  #1  and  A  #2,  Class  I,  Method  1,  61/ 
frequency 
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a=0. 162567 
b=0. 190477 


A 

a  a 


f=0.205908 

g=0.198013 


c=0.181908  a  a  .  h=0.215354 

d=0. 187383  A  i=0.216628 

e=0. 202820  deed 

/frvfc-X 

e  g  f  f  g  e 

e  f  hi  i  h  f  e 

A  9_)f 1 A '  A  9-X 

d  e  f  h  i  i  h  f  e  d 

A  cA f  A  h-)f '  A  c;X 

abcegf  f  geeba 

L,  -V-e  -V-e  -V-d  -V-a  -A 


Figure  4.10:  Icosahedron  chord  lengths,  R  =  1 


Figure  4.11:  Sub-icosa  (panel)  types 
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(K.LM1.1) 

(N.MH2.1) 


2 


Figure  4.16:  Icosa  (1, 1)  showing  subicosa  location  (2, 1) 


Sublcosa  Triangle  in  Local  (>^M.yNM)  Coordinates 


Figure  4.17:  Subicosa  triangle  (2, 1)  in  local  (xNm,Unm)  coordinates 
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Figure  4.18:  Icosa  (1, 1)  showing  subicosa  location  (2, 2) 


Figure  4.21;  Subicosa  triangle  (2, 


(K.LM1.1) 

(N,M)=(3,1) 


Figure  4.22:  Icosa  (1, 1)  showing  subicosa  location  (3, 1) 


Sublcosa  Triangle  in  Local  (^,M.yNM)  Coordinates 


Figure  4.23:  Subicosa  triangle  (3,1)  in  local  ( xnMjUnm )  coordinates 
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Figure  4.26:  Icosa  (1, 1)  showing  subicosa  location  (3,3) 


^NM  y  NM 

Figure  4.30:  Subicosa  and  equilateral  triangles 


a=0. 162567 


f=0.205908 


b=0. 190477  A  g=0. 198013 

c=0.181908  a7  a  h=0.215354 

d=0. 187383  i=0.216628 

e=0.202820  deed 

/f  e"X  . 

e  g  f  f  g  e 

Ar  9")f 1  ~¥c 1  9_X 

d  e  f  hi  i  h  f  e  d 

k  c_)f f  ~¥r  h_)f f  C_X 

abceg  ffgecba 

Z-a-V-d-V-e-^-e-V-d-V-aA 


Figure  4.31:  Icosahedron  cord  lengths  indicating  base-lengths 
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Figure  4.33:  Equilateral  triangle  pertaining  to  determination  of  the  rotation  angle  <pnm 


Figure  4.34:  Local  and  reference  coordinates  pertaining  to  determination  of  the  rotation 
matrix 
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Table  4.3:  Rotation  angles  ipuM  for  icosa  ( K ,  L)  =  (1, 1) 


Rotation  angle^for  Icosa  (K,L)=(2»1) 

1 

2 

3 

4 

5 

6 

7 

a 

9 

10 

11 

180,0000 

181.0433 

178.9567 

* 

-80.2441 

119.2178 

180.0000 

240.7822 

60,2441 

176.8323 

59.0210 

0.0000 

-59.0210 

121.5295 

183.3674 

-66.6925 

-3.4444 

-83,1128 

118.6039 

180.0000 

241.3961 

63.1128 

■ — . ; 

— 

-58.3393 

124.6141 

69.3932 

Table  4.4:  Rotation  angles  for  icosa  (K,  L )  =  (2, 1) 
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Rotation  angle  for  icosa  (K,L] 

=(2,2) 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

0.0000 

1.0433 

180.0000 

-1.0433 

119.7559 

-60.7822 

0.0000 

60.7822 

240.2441 

-3.3674 

58.4705 

239.0210 

180.0000 

120.9790 

-58.4705 

3.3674 

113.3075 

176.5556 

_ ! 

116.8872 

-61.3961 

0.0000 

183.4444 

246.6925 

pwi 

110.6068 

249.3932 

Table  4.5:  Rotation  angles  <pnm  for  icosa  (K,  L)  =  (2,2) 


Rotation  angle^^or  icosa  (K,L)=(3,1)  | 

N\M 

1 

2 

3 

4 

5 

6 

7 

8  , 

9 

10 

11 

1 

0.0000 

2 

21.9580 

180.0000 

-21.9580 

: 

3 

144.9632 

-45.7860 

0.0000 

45.7860 

215.0368 

4 

23.8554 

78.9418 

248.9854 

180.0000 

111.0146 

-78.9418 

-23.8554 

5 

142.4488 

200.7306 

132.7726 

-52.7103 

0.0000 

52.7103 

227.2274 

159.2694 

217.5512 

6 

142.0906 

-39.7757 

135.2281 

70.3892 

245.5333 

180.0000 

114.4667 

-70.3892 

224.7719 

39.7757 

217.9094 

Table  4.6:  Rotation  angles  for  icosa  (K,  L )  =  (3, 1) 
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(K,L)=(1,1) 

(K,L)=(2,1) 


Figure  4.35:  Four  Basic  Icosas 
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Figure  4.38:  Triangle  Type  1 


Figure  4.39:  Triangle  Type  2 


Figure  4.40:  Triangle  Type  3a 


Figure  4.41:  Triangle  Type  3b 


I 

I 


i 

h  i  ar 

i  £ 


Figure  4.42:  Triangle  Type  4 


Figure  4.43:  Triangle  Type  5a 
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Figure  4.45:  Triangle  Type  6 
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Single  fcosa  -  Top  View 

(N.MH1.1) 


Figure  4.46:  Subicosa  (1?  1, 1?  1)  -  top  view 


Sublcosa  and  Equilateral  Triangles  In  Reference  Cbordfnates 


Figure  4.47:  Subicosa  (1,  1, 1)  and  respective  equilateral  triangle  in  reference  coordinate 
system  (x'NM,y'NM) 
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Single  Icosa  -  Top  View 
(K.LM1.1) 

(N,M)={2,2) 


Figure  4.48:  Subicosa  (1, 1, 2, 2)  -  top  view 


Sublcosa  and  Equilateral  Triangles  in  Reference  OfgM.y*NM)  Coordinates 


Figure  4.49:  Subicosa  (1}1, 2, 2)  and  respective  equilateral  triangle  in  reference  coordinate 
system  (x’NM,y'NM) 
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Single  Icosa  -  Top  View 

(K.tW.1) 


Figure  4.50:  Subicosa  (1*  1,  25 1)  -  top  view 


Sublcosa  and  Equilateral  Triangles  in  Reference  Coordinates 


Figure  4.51:  Subicosa  (1,  Ij2?  1)  and  respective  equilateral  triangle  in  reference  coordinate 
system  {x'NM)y'NM) 
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Single  Icosa  -  Top  View 

(K.lhi.D 

(N,M)=(2,3) 


Figure  4.52:  Subieosa  (1, 1,2,3)  -  top  view 


Sublcosa  and  Equilateral  Triangles  in  Reference  (^.y^)  Coordinates 


Figure  4.53:  Subieosa  (1,1, 2,3)  and  respective  equilateral  triangle  in  reference  coordinate 
system  (x'NM,y'NM) 
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Single  Icosa  -  Top  View 
0CLM1.1) 

{N,M}={3,3) 


Figure  4.54:  Subicosa  (1?  13  3*  3)  -  top  view 


Sublcosa  and  Equilateral  Triangles  in  Reference  (>^M,y^M)'Coordinates 


Figure  4.55:  Subicosa  (1, 1, 3, 3)  and  respective  equilateral  triangle  in  reference  coordinate 
system  {x'NM,y'NM) 
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ive  equilateral  triangle  in  reference  coordinate 


07 


Single  Icosa  -  Top  View 

(K.L)=(1,1) 

(N,M}={4,4) 


Figure  4.60:  Subicosa  (1, 1,4,4)  -  top  view 


Sublcosa  and  Equilateral  Triangles  in  Reference  Coordinates 


Figure  4.61:  Subicosa  (1,1, 4, 4)  and  respective  equilateral  triangle  in  reference  coordinate 
system  (x'NM,y'NM) 
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Figure  4.64:  Equilateral  triangle  pertaining  to  location  of  array  elements  in  reference  coor¬ 
dinate  system 
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m=1,2,...1(nmMax+1-n)  dy =3,022" 
(nIm)=(1l1) - ►  ~*\  I"*-  | 


Figure  4.65:  Array  elements  in  reference  coordinate  system  (x nm^nm) 


Panel  Type 

(N,M) 

Equilateral  triangle  side  (m) 

No  of  Elements 

R=tm 

R=5m 

R=5m 

i 

(1.D 

0,1521 

0.7607 

8 

36 

2 

(2.2) 

0.179 

0,8948 

10 

55 

3 

(2.1) 

0,1781 

0,8906 

10 

55 

4 

(3.3) 

0.2027 

1.0133 

12 

78 

5 

(3,1) 

0,198 

0.9901 

11 

66 

6 

(4,4) 

0.2154 

1,0768 

13 

91 

Table  4,7:  Number  of  rows  and  elements  for  6  panel  types 
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Equilateral  Subarray 


Figure  4.66:  Equilateral  subarray  in  reference  (x'NM,y'N M)  coordinates,  R  =  5m,  Type  1 


Equilateral  Subarray 


-0.6  -0.4  -0.2  0  0.2  0.4  0.6 


Wm> 


Figure  4.67:  Equilateral  subarray  in  reference  (x'NM,  y'NM)  coordinates,  R  =  5m,  Type  2 


113 


Equilateral  Subarray 


Chapter  5 

Geodesic  Dome  Array  Beamsteering 
and  Radiation  Patterns 


5.1  Element  pattern 


5.1.1  Element  pattern  in  reference  coordinates 


The  element  pattern  is  originally  given  in  reference  coordinate  system  (x'NM,  y'NM,  z'NM)  (see 
Figure  4.65),  which  for  CP  field  can  be  written  as 


gnm(Qnm’  4>nm)  —  E0yJ cos  6'NM  e±^NM— f=  (6NM  ±  (5-1) 

where  ±  refers  to  the  LHCP  and  RHCP  wave  fields.  Here,  (0'NM,  4>'nm)  8X6  the  observation 
angles  in  reference  coordinates.  Since  every  element,  (nm),  on  a  given  subicosa  (subar¬ 
ray/facet/panel)  sees  the  far  field  observation  point  at  the  same  angle,  in  (5.1),  we  replaced 
55 element”  subscripts  (nm)  by  ”panel”  subscripts  (NM),  Furthermore;  in  (5.1)?  (Omm^nm) 
are  respective  unit  vectors  as  shown  in  Figure  5.1,  and 


E0 


I  Co  /47Tj4ceii 

V  4tt  V  A2 


(5.2) 


5.1.2  Element  pattern  in  local  coordinates 

Next,  we  transform  this  element  field  pattern  into  local  coordinates  (x^m,  Vnm,  znm)  by 
using  following  identities: 


®NM 

‘  1  0  ' 

&NM 

_  Qnm  . 

.  4>nm  . 

x'nm 
Vnm 
.  znm  . 


COS  if  MM  f  MM  0 

XtfM 

— 

—  sin  pmm  cos  pmm  0 

Vnm 

1 

o 

o 

h-» 

i _ 

znm 

(5.3) 

(5.4) 
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&NM  ~  @NM  (5-5) 

<I>'nm  ~  &NM  ~  <PNM-  (5-6) 

Thus,  in  local  coordinates  (see  Figure  5.2),  the  element  patterns  becomes 

- -  1 

&nm(9nm,  4>nm)  —  Eo\f  cos  9nm  e±j<t>NM— ^  {Qnm  ±  J&nm)-  (5-7) 

5.1.3  Element  pattern  in  global  coordinates 

The  respective  local  unit  vectors  (x-nMiYnm,  zatm)  are  related  to  the  unit  vectors  of  the 
global  coordinate  system  (x,  y,  z)  by  the  following  rotation  matrix  equation: 

Xjvm  _  x 

y  nm  =  A-nm  y  (5-8) 


where  A  is  a  rotation  matrix  and  is  given  by 

cos  a nm  cos  Pnm  cos  a nm  sin  P nm  —  sin  ocnm 
ANm  =  sin  (3 nm  cos  Pnm  0  .  (5.9) 

sin  ocnm  cos/3nm  sin  &nm  sin  Pnm  cos  q^nm 

From  Appendix  B.3,  the  following  two  relations  are  obtained: 

cos  9nm  —  sin  c»-nm  sin  9  cos(cf)  —  Pnm)  +  coso^vf  cos  9- . . ,-(5.10) 

cos  aNM  sin  9  cos {p  —  Pnm)  —  sin  cunm  cos  9  .  . 

cot  <pNM  —  - ■  .  ■  , - n - v - •  (5.11) 

sm9sm(4>  -  Pnm) 

Similarly  the  local  unit  vectors  {Qnm,  4>nm)  can  be  transformed  into  global  coordinates 


where 


O.NM  = 


~aNM  £>nm  Q 

—1>NM  —O-NM  J  I  (j) 


cos  9  sin  aNM  cos{<p  —  Pnm)  —  sin  9  cos  a atm 


(5.12) 


(5.13) 


t>NM  = 


sin  aNM  sin(0  —  Pnm) 


(5-14) 


Consequently,  in  (5.7),  as  in  Appendix  B.6  for  an  array  element,  the  factor  {Qnm±]4>nm) 
becomes: 

Qnm  ±  J4>nm  =  ~{Q  ±  ]4>)  «±jVnm  (5-15) 


where 


tan  vnm  = 


bm£  _ _ sin  qjvm  sin(<ft  —  Pnm) _ 

clnm  cos  9  sin  onm  cos(o  —  Pnm)  ~  sin  9  cos  a  nm 


(5.16) 
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Now  we  can  rewrite  the  element  pattern  relation  (5.7)  in  global  coordinates  as 

enm(9,  <f)  =  —Eq\J cos  9Nm  (9  ±  j$)  e±jUNM  e±j*™  (5.17) 


where  cos  O^m  and  4>nm  are  expressed  in  terms  of  (0,  <fr)  by  (5.10)  and  (5.11).  Here  the  phase 
ipnm  accounts  for  the  difference  in  propagation  path  lengths  from  the  radiating  element  to 
the  far  field  point  and  the  global  origin  to  the  far  field  point.  In  reference  to  Figure  5.2,  ipnm 
is 

$ nm  =  kKnm-t .  (5.18) 

The  radial  vector  specifying  the  location  of  the  array  elements  in  global  coordinates,  Fig¬ 
ure  4.32  and  Figure  4.64,  is 

Rnm  ~  RiW  d"  ^nm  (5.19) 

and  is  defined  in  (4.78)  as 

Rnm  =  xxnm  +  yyn  m  zznm  (5.20) 

where  (xnm,ynm,znTn)  are  the  coordinates  of  the  array  elements  in  global  coordinates.  Since, 

r  =  xsin5cos0-t-ysin0sin0  +  zcos0  .  (5.21) 


from  (5.18)  using  (5.19)  and  (5.21)  we  get 

ip  =  k(xnm  sin  6  cos  <f>  +  ynm  sin  6sin<p  +  znm  cos  6) .  (5.22) 

Notice  that  in  (5.17)  we  re-introduced  subscripts  (n,  m),  i.e.,  we  replaced  by  enm  because 
factor  ipnm  depends  on  (ram). 


5.2  Far-field 

Array  far-field  (in  global  coordinates): 

pi{ujt-kr) 

E(r;  t )  =  E(0,  </>) - —  (5.23) 

r 

where 

E(0,  <t>)  =  6Ee(0,  <t>)  +  4>E^(0,  <f>)  =  £Enm(0,  <fi)  (5.24) 

nm 

and,  the  field  due  to  a  single  element  in  the  array  environment  is 

Enm($>0)  =  finm  Cnmifi 0>  <fe)  ®nm($i^)-  .  (5.25) 

with  enm  given  by  (5.17). 
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5.2.1  On-Off  switch: 

In  (5.25),  the  Kronecker  delta  5nm  defines  the  active  sector  on  the  array.  In  other  words  it 
is  a  switch  turning  on  and  off  the  array  elements  as  beam  scans.  We  used  two  criteria  to 
determine  the  active  sector:  (a)  element  level  and  (b)  panel  or  subarray  level. 

(a)  element  level:  In  this  case  all  elements  within  the  cone  angle  amax  are  on,  while 
all  others  are  turned  off,  Figure  5.3  and  Figure  5.4.  The  axis  of  the  cone  coincides  with  the 
main  beam  direction.  This  criteria  gives  smooth  edges  to  the  active  ’’spherical”  sector  which 
means  it  can  potentially  produce  a  low  sidelobe  radiation  pattern  if  required.  The  drawback 
of  this  criteria  is,  however,  that  since  some  of  the  elements  located  on  the  ’’edge”  panels 
(subarrays)  will  be  turned  on  and  some  will  be  turned  off,  there  will  be  degradation  of  the 
S/N  or  G/T  on  the  panel  level  in  the  receive  case  or  reduced  efficiency  in  the  transmit  case. 
We  assumed  that  each  panel  has  one  output  per  beam  and  that  all  subarray  elements  are 
fed  by  the  corporate  RF  manifold.  In  addition,  this  scheme  requires  more  computational 
power  in  the  beamsteering  algorithm  than  the  scheme  (b)  which  for  a  large  array  may  be 
undesirable.  It  may  be  interesting  to  point  out  that  the  edge  elements  on  different  panels 
see  the  far  field  at  the  peak  of  the  beam  at  slightly  different  angles  (within  few  degrees) 
with  respect  to  the  local  normal.  For  example,  if  we  set  amax  —  60°,  the  local  observation 
angles  of  the  edge  elements  will  have  different  values  in  the  vicinity  of  60°  plus/minus  a 
couple  of  degrees  so  the  edge  circle  may  exhibit  a  small  discontinuity  as  it  crosses  from  one 
to  neighboring  panel.  Thus,  we  define 

r  _  J  1  if  7nm  ^  amax  /c  r>c\ 

°nm  \  0  otherwise  (5>26) 

where  'ynm  is  the  angle  between  the  axis  of  the  cone  and  the  radial  line  passing  through  the 
element  ( nm )  as  shown  in  Figure  5.3.  To  determine  7 nm  we  write 


COS  'Ynm  —  f0  ■  Tnm  (5.27) 

Here,  the  unit  vector  in  direction  of  (n,  m)  element  is 


ffnn 

Rnn 


(5.28) 


where  Rnm  is  defined  by  (5.20)  and  the  distance  between  the  element  (nm)  and  the  origin 
of  the  global  coordinate  system  is 

Rnm  =  Jxln  +  V L  +  4n-  (5-29) 

Furthermore  in  (5.27),  fo  is  the  unit  vector  in  the  main  beam  direction  (0o,  (f>o),  thus 

r  o  =  x  sin  0O  cos  </>0  +  y  sin  0O  sin  <f>0  +  z  cos  90 .  (5.30) 

Substituting  (5.20)  and  (5.29)  into  (5.28)  and  then  (5.28)  and  (5.30)  into  (5.27),  we  get 

cos  7„m  =  sin  0O  cos  </>0  +  sin  0O  sin  <£0  4-  cos  0O .  (5.31) 

ftnm  *Wim 
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(a)  panel  level:  In  this  case  Snm  in  (5.25)  is  replaced  by  SNM  where  (NM)  refers  to  the 
subarray  containing  element  ( nm ).  In  this  scheme  the  active  sector  on  the  dome  consists  of 
a  fixed  number  of  active  panels  while  all  other  panels  are  completely  off.  As  a  beam  scans 

across  the  sky,  some  of  the  edge  panels  will  be  switched  off  while  some  panels  that  were 

off  will  be  switched  on.  It  is  clear  that  this  scheme  is  simpler  than  scheme  (a),  therefore 
it  require  less  computational  power  for  the  beamsteering.  Also,  compared  to  scheme  (a), 
there  is  no  additional  RF  loss  in  the  edge  panels  because  each  individual  panel  has  either  all 
elements  turned  on  or  off.  The  drawback  of  this  scheme  is  that  the  active  sector’s  edges  are 
rough  and  because  of  that  the  beam  shape  won’t  be  uniform.  Also  it  may  not  give  very  low 
sidelobes  especially  for  small  active  areas  on  the  dome.  In  this  case, 

x  _  J  i  ^  1 nm  <  ^■TTXCLX  /r 

=  \  0  otherwise  (5'32> 

where  jnm  is  the  angle  between  the  beam  axis  and  the  radial  line  passing  thru  the  origin  of 
the  local  coordinate  system  Onm(Rnm,  &nm,  Pnm)-  Thus  all  panels  with  Onm  falling  inside 
the  cone  angle  amax  will  be  turned  on  while  all  others  will  be  off.  In  reference  to  Figure  5.5 
and  Figure  5.6,  7 nm  is 

cos  7^  =  r0  •  tnm  (5.33) 

where 

Tnm  —  x  sin  oinm  cos  Pnm  4-  y  sin  ocnm  sin  Pnm  +  z  cos  oinm-  (5.34) 

Substituting  (5.30)  and  (5.34)  into  (5.33),  we  get 

cos  7 nm  —  sin  6q  cos  <fio  sin  cx-nm  cos  Pnm  +  sin  60  sin  <fio  sin  ocnm  sin  Pnm  +  cos  8q  cos  ocnm- 

(5.35) 

In  (5.35),  see  Figure  5.5,  amax  determines  the  size  of  the  active  sector.  Typically,  amax  <  60°. 
In  this  case  the  array  elements  on  the  edge  panels  will  see  the  beam  observation  point  at 
local  angles  60°  or  less,  but  not  more  than  amax.  Local  observation  angles  are  equivalent  to 
planar  subarray  scan  angles  measured  from  local  broadside.  As  is  known,  the  planar  array 
impedance  and  the  radiation  and  polarization  properties  rapidly  degrade  after  60°  scan,  so 
it  is  not  advisable  to  set  amax  greater  than  60°.  So,  in  this  scheme,  the  maximum  scan  angle 
of  the  subarray  is  equal  or  less  than  amax.  This  is  an  important  parameter  to  know  in  the 
design  procedure  of  the  subaxray. 

The  panel  level  scheme  also  has  the  advantage  that  the  numbers  of  beams  provided  by 
subarray /panel  can  be  doubled  or  quadrupled  with  very  small  additional  cost.  To  elaborate 
on  this,  we  assume  that  each  panel’s  BFN  consists  of  two  transmit  and  two  receive  RF 
manifolds.  This  means  that  the  panel  can  produce  four  simultaneous  independent  beams 
(two  transmit  and  two  receive  beams).  However,  not  all.  satellite  links  require  frill  size 
aperture.  Satellites  in  deep  space  GEO  or  HEO  may  need  full  dome  aperture,  but  MEO 
and  LEO  satellites  may  need  only  a  fraction  of  the  full  size  aperture  (six  or  twelve  panels 
for  example).  The  geodesic  dome  antenna  can  provide  ’’gain  on  demand”,  i.e.,  the  EIRP 
for  transmit,  beams  and  G/T  for  receive  beams  can  be  adjusted  to  meet  specific  satellite  link 
requirements.  This  is  done  by  varying  size  of  the  active  sector  as  illustrated  in  Figure  5.9 
for  two  independent  transmit  beams,  TX1-A  and  TX2-A. 
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The  geodesic  dome  array  based  on  icosahedron  structure  with  Qu  subdivision  frequency 
consists  of  675  subarrays/panels.  As  already  mentioned,  panels  for  a  given  beam  in  an  active 
sector  are  turned  on  while  all  others  are  off.  Panel  beam-ports  that  are  off  can  be  used  to 
form  an  additional  four  or  more  beams.  This  can  be  done  simply  by  adding  four  single-pol 
double-throw  RF  switches  at  the  panel  output  beam  ports  and  four  1:675  RF  manifolds  as 
illustrated  in  Figure  5.8  for  the  special  case  where  we  double  the  number  of  transmit  and 
receive  beams.  The  two  additional  transmit  beams  are  designated  by  TX1-B  and  TX2-B  in 
Figure  5.9.  It  is  important  to  mention  that  the  active  areas  of  two  beams  from  the  same 
subarray  beamport  can  not  overlap.  The  TX1-B  and  TX2-B  are  completely  independent 
beams,  however  the  active  sector  of  TX1-B  can  not  overlap  with  TX1-A  and  active  sector 
of  TX2-B  can  not  overlap  with  TX2-A.  Thus,  in  this  case  the  geodesic  dome  can  support 
four  transmit  beams  and  four  receive  beams,  two  transmit  and  two  receive  beams  being 
completely  independent  while  the  other  four  beams  (two  Tx  and  two  Rx)  are  conditionally 
independent  . 

In  practice,  since  the  satellites  are  roughly  uniformly  distributed  in  a  hemisphere,  using 
the  antenna  resource  management  computer  (smart  time  and  space  satellite  link  management 
and  scheduling)  the  conditional  independence  can  be  met  most  of  the  time.  In  the  case  this 
condition  can’t  be  satisfied  at  first,  beam  re-assignment  on  the  fly  may  be  required  (switching 
with  another  beam/satellite).  If  the  independence  condition  between  the  two  beams  is  still 
not  met,  then  one  can  resolve  this  situation  by  moving  one  or  both  relevant  active  sectors  off 
conical  axis  (away  from  each  other).  Consequently,  the  two  beams  in  conflict  are  no  longer 
broadside  but  are  skewed  with  respect  to  their  cone  axis  (broadside)  and  therefore  will  suffer 
small  scan  losses.  The  antenna  resource  management  computer  can  compensate  this  losses 
by  adding  one  or  two  panels  to  their  active  areas  needed  at  broadside. 


5.2.2  Element  voltage  excitation  coefficients 

In  (5.25),  cnm  are  the  element  voltage  excitation  coefficients 

Cnm(do,  to)  =|  cnm(6o,  to)  I  (5.36) 

where  |  cnm  |  is  the  desired  aperture  illumination  function. 

From  (5.25)  and  (5.17)  it  is  seen  that  to  point  the  beam  in  (0O,  to)  direction,  Figure  5.7, 
Xnm  must  be 

Xnm{6o,  to)  =  -[Vw(0O,  to)  ±  VNm(8 0.  to)  ±  tnm{Q 0,  to)]-  (5.37) 


5.3  Array  Gain 

The  array  gain  is 


G(e,t)  = 


E(r)  •  E*(r) 

Co 


(5.38) 


where  free  space  impedance  Co  =  376.7  ohms,  and  Pq  is  isotropic  source  power  density,  i.e., 


Bo  = 


Enm  ti 


Cnm 


Airr2 


(5.39) 
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Thus, 


G{e,<t>)  = 


4^  |  Enm  $nmCnm&nm{& j  $) 


Co 


Enm  ^nm 


Cnm 


From  here,  peak  array  gain  is 


G(0 0,to)  = 


4tt[E. 


nm  4m 


£nm 


i(Go,<Po)  \f 


Co 


EZnm  4m 


^nm 


(5.40) 


(5.41) 


It  is  important  to  keep  in  mind  that  in  (5.41),  for  example  for  icosahedron  with  subdivision 
frequency  u  =  6,  using  Class  I,  Method  1,  subdivision  scheme 

3  5(10)  6  2N— 1  nmMax  nmM ax+l—n 

£=££££  £  £  (5.42) 

nm  K= 1  i=l  AT=1  M=1  n=  1  m=l 

where  as  already  discussed  (/CL)  denotes  icosa  of  the  icosahedron,  (NM)  denotes  the  sub- 
icosa  or  subaxray  panel  (facet)  and  (nm)  denotes  the  radiating  element  on  this  subarray.  If 
K= 1  or  2, 1=1  to  5,  and  if  K—2  then  L  varies  from  1  to  10. 
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Figure  5.1:  Observation  angles  {8'mmA'nm)  in  reference  coordinates 


Figure  5.5:  Active  array  sector  -  panel  level 


Figure  5.6:  Geometry  pertaining  to  evaluation  of  jmm 


TX2-B 


TX1-B 


Figure  5.9:  Geodesic  dome  showing  active  sectors  of  four  simultaneous  transmit  beams  (four 
receive  beams  not  shown) 
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Chapter  6 

Computer  Codes 


6.1  Code  1 

6.1.1  Input 

Based  on  the  theory  presented  above,  we  have  developed  a  MATLAB  code  ”boris5.m”  that 
evaluates  all  relevant  geometrical  parameters  necessary  for  the  construction  geometry  and 
beamsteering  of  the  geodesic  dome  antenna.  For  a  given  dome  radius,  this  code  needs  to  be 
run  only  once.  The  code  saves  all  output  data  into  file  ”dataGD5.mat”. 

Input  parameters  can  be  entered  via  the  GUI  shown  in  Figure  6.1.  For  a  selected  icosa 
(K,L),  and  subicosa  ( N,M ),  and  radiating  element  ( n,m )  the  code  produces  a  number 
of  plots  that  may  be  useful  for  the  design  of  a  geodesic  dome.  For  example  it  plots  the 
icosa  and  subicosa,  and  the  respective  radiating  elements  on  a  3D  dome.  By  zooming  into 
a  subicosa  one  may  see  all  radiating  elements  where  the  selected  element  (n,  m)  is  fleshed 
out.  Furthermore,  the  selected  subicosa  is  plotted  in  the  local  coordinate  system  as  well  as 
in  the  reference  coordinates.  The  respective  triangle  coordinates  are  printed.  Similarly,  the 
equilateral  triangles  and  array  elements  are  plotted  in  reference  coordinates. 

There  is  one  single  input  parameter,  i.e.,  array  radius  R.  The  code  is  valid  for  subdivision 
frequency  6,  class  I.  The  code  saves  output  data  if  ’’Save  Data”  button  is  checked  on. 

6.1.2  Output 

The  output  data  file  saves  the  following  parameters: 

XI  Y1  Z1  X2  Y2  Z2  X3  Y3  Z3  =(x,y,z)  coordinates  of  subicosa’s  (GD  panels).  It 
evaluates  equation  (4.1).  The  4D  arrays  are  of  the  size  DIM(K,L,N,M) 

XO  YO  ZO  =(xq,  ?/o,  zq)  coordinates  of  the  origin  of  local  coordinate  system  (xaw,  Unm,  znm)- 
See  equation  (4.15).  DIM(K,L,N,M) 

R0  ALPHA  BETA=(R,ym,  &nm,@nm)  coordinates  of  the  origin  of  local  coordinate  sys¬ 
tem  (xjvm, VnMi znm)-  Eqs.  (4.16)  to  (4.18).  DIM(K,L,N,M). 

LX  MX  NX  LY  MY  NY  LZ  MZ  NZ  =directional  cosines  of  local  (xjvm,  Vnm,  znm)  axes, 
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DIM(K,L,N,M).  Eqs.  (4.19),  (4.21),  (4.23)  and  (4.24). 

TYPE  =  panel  (subicosa)  type,  DIM(N,M),  see  Figure  4.11.  . .  .  .  _ 

PhiP  =rotation  angle  <pnm  in  radians,  (angle  between  and  x'NM),  DIM(K,L,N,M). 
Equation  (4.32). 

aE  =  equilateral  triangle  side  length,  DIM(K,L,N,M),  See  Section  4.3.2. 

XCEm  YCEm  =  center  of  mass  for  the  equilateral  triangle  in  reference  coordinate  system. 
Wnm^nm^nm)^  DIM(K,L,N,M),  see  Section  4.3.2.. 

nmMax  =Number  of  array  rows  in  a  panel,  DIM(N,M),  equation  (4.63) 

xEL  yEL  zEL  =  array  element  location  in  global  (x,  y,  z )  coordinate  system, 
DIM(K,L,N,M,n,m) ,  eq.  (4.74). 

The  user  friendly  code  was  generated  and  tested  in  Matlab  6.5,  Release  13.  To  run  the 
code  the  user  should  go  to  the  command  window  and  type  boris5.m  -  Enter.  Then  the  GUI 
window  will  pop  up.  Set  the  input  parameters  as  desired  and  press  the  run  button.  It  is 
suggested  that  the  user  should  first  run  the  code  with  default  input  parameters  and  compare 
the  output  data  with  the  original  values  and  plots  in  Figure  6.2  to  Figure  6.9. 

By  selecting  different  (K,  L,  N,  M,n,m)  one  may  observe  specific  panel  and  element 
locations  on  the  dome.  Coordinates  and  all  relevant  geometrical  parameters  can  be  read 
either  from  the  plots  for  selected  ( K,L,N,M,n,m ),  or  by  going  to  the  Matlab  workspace 
to  read  the  above  output  parameters  which  are  written  in  a  multidimensional  array  form. 
Leave  the  ’’Save  Data”  button  off  all  the  time  except  when  the  output  data  must  be  stored 
in  output  file  dataGD5.mat.  In  principle,  for  given  sphere  radius  R,  the  code  needs  to  be 
executed  only  once. 

6.2  Code  2 

The  second  Matlab  code  ”Tomasic5.m”  evaluates  the  required  element  phasing  for  a  beam 
in  the  desired  (0o,0o)  direction;  radiation  patterns  are' also  generated.  The  code  reads  the 
input  data  file  ”dataGD5.mat”  that  contains  all  necessary  array  geometry  for  beamsteering 
and  calculation  of  the  far-field.  The  code  plots  a  full  dome  array  with  all  radiating  elements. 
The  elements  in  an  active  sector  are  plotted  in  red  while  elements  turned  off  sire  in  green. 
Using  the  Matlab  zoom  command  and  object  rotation,  the  user  can  interactively  focus  on 
specific  area  of  interest  around  the  dome.  Since  there  are  tens  of  thousands  of  elements  in 
the  array,  the  interactive  plotting  routine  is  unfortunately  slow,  but  the  graphics  are  useful 
for  understanding  and  visualizing  the  geodesic  dome  array  in  action. 

6.2.1  Input 

Similarly  as  in  Code  1,  the  input  parameters  can  be  entered  via  the  GUI  shown  in  Figure  6.10. 
The  first  column  in  Figure  6.10  specifies  the  selected  panel  and  element  on  the  dome.  The 
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selected  panel  and  element  will  be  fleshed  out  in  red  in  the  dome  plot.  The  second  column 
in  Figure  6.10  sets  the  array  radius,  frequency,  the  beam  pointing  direction  ( do ,  4>a),  the 
active  sector  angle  amax,  and  the  ^=constant  plane  in  which  the  radiation  pattern  will  be 
evaluated  and  plotted.  If  the  user  wants  to  capture  the  peak  of  the  beam,  then  one  must 
set  <f)  =  4>q.  Finally,  third  column  in  GUI  Figure  6.10,  defines  two  aperture  illuminations, 
uniform  and  s/cosd.  If  some  other  low-sidelobe  aperture  taper  is  needed  then  It  is  easy  to 
modify  the  program  accordingly.  Then  we  have  a  choice  of  two  criteria  for  the  active  sector: 
criteria  based  on  the  element  or  on  the  panel  level.  Furthermore,  we  have  two  choices  for  the 
polarization,  circular  or  elliptical;  and  a  of  choice  whether  we  want  to  plot  both  the  radiation 
pattern  and  the  full  geodesic  dome  or  just  the  pattern. 

6.2.2  Output 

The  GUI  in  Figure  6.10  shows  the  default  input  parameters.  Figure  6.11  shows  the  respective 
active  sector  on  the  dome.  Figure  6.11  and  Figure  6.12  zoom  in  on  the  selected  panel,  while 
Figure  6.14  plots  the  radiation  pattern  for  input  of  Figure  6.10. 

Figure  6.15  shows  the  same  GUI  except  the  active  sector  is  defined  on  the  panel  level. 
For  an  active  sector  with  amax  —  60°,  and  (00, 4>o)  =  (0, 0),  the  sector  is  in  a  pentagon  shape 
as  shown  in  Figure  6.16  and  in  zoom  in  version  in  Figure  6.17.  The  respective  radiation 
pattern  is  shown  in  Figure  6.18.  Figure  6.20  and  Figure  6.21  are  output  plots  with  input 
parameters  given  in  GUI  on  Figure  6.19.  This  case  is  the  same  as  that  pf  Figure  6.15 
except  here  Qq  =  45°.  Comparing  Figure  6.18  and  Figure  6.21  we  see  that  the  radiation 
pattern  is  practically  independent  of  scan,  which  is  an  expected  and  desirable  feature  of  the 
spherical  arrays.  Figure  6.22  to  Figure  6.25  shows  the  same  except  in  this  case  we  apply 
panel  level  active  sector  criteria.  The  active  sector  for  this  case  is  shown  in  Figure  6.23  with 
the  respective  radiation  pattern  shown  in  Figure  6.25.  Notice,  by  comparing  Figure  6.21  and 
Figure  6.25,  we  see  that  for  a  large  active  sector  there  is  not  much  difference  in  the  radiation 
patterns  between  the  two  criteria. 

For  LEO  satellites,  the  EIRP  and  G/T  link  requirements  are  much  lower  than  for  GEO 
satellites,  and  therefore  the  active  sector  on  the  dome  could  be  much  smaller  than  in  previous 
examples.  For  example,  the  active  sector  of  an  equivalent  4m  dish  is  shown  in  Figure  6.26. 
This  spot  will  give  roughly  G/T  =  11  dB/K  enough  for  most  low-flying  birds.  The  active 
spot  has  a  smooth  contour  which  indicates  that  we  used  element  level  criteria  to  determine 
the  active  sector.  The  respective  radiation  pattern  is  shown  in  Figure  6.27.  Figure  6.28  to 
Figure  6.35  show  the  same  in  scan  intervals  in  elevation  of  5  deg.  Note  the  radiation  pattern 
and  peak  gain  as  beam  sweeps  from  6q  =  37.38  deg  to  57.38  deg  are  fairly  scan  independent. 

To  see  the  rough  edge  effects  on  radiation  pattern  and  peak  gain  we  run  the  same  case 
but  the  active  sector  is  determined  based  on  whole  panel  criteria.  As  seen  from  Figure  6.36 
to  Figure  6.45  the  peak  gain  is  fairly  constant  while  sidelobes  are  little  higher  than  in  the 
previous  smooth-edge  case. 
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Figure  6.1:  GUI  for  input  parameters 
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Icosahedral  Dome 


(k,lhi.i) 

(N,M)=(1,1) 


Figure  6.2:  Icosahedral  dome  in  perspective  indicating  selected  icosa 


Figure  6.3:  Zoom  in  of  Figure  6.2 
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Single  Icosa 

(K.L)=(1,1)  ■  z 

(N,M)=(1,1) 


Figure  6.5:  Single  Icosa 
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Figure  6.7:  Equilateral  subarray  indicating  (n,m)=(l,l 


Sublcosa  and  Equilateral  Triangles  in  Reference  (^M,y^M)  Coordinates 


Figure  6.8:  Subicosa  ( K ,  L,  N,  M)  —  (1, 1, 1, 1)  and  respective  equilateral  triangle  in  reference 
coordinates 


Sublcosa  Triangle  in  Local  0^M.yNM)  Coordinates 


Figure  6.9:  Subicosa  ( K ,  L,  N ,  M)  —  (1, 1, 1, 1)  in  local  coordinates 
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Figure  6.10:  GUI  Code  2 


Figure  6.11:  Full  dome  showing  amax 


60°  active  sector  in  red 


Panel  Radiation  Pattern 


Figure  6.14:  Radiation  pattern  pertaining  to  input  parameters  of  Figure  6.10  -  default  input 
values  - .  r 
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Figure  6.17:  Zoom  in  of  Figure  6.16 


Figure  6,18:  Radiation  pattern  pertaining  to  GUI  of  Figure  6.15 
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Panel  Radiation  Pattern 


Figure  6.21:  Rad  pattern,  45  deg  scan 
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Figure  6.22:  GUI  45  degree  scan  -  active  sector  ’’panel” 
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Figure  6.23:  Active  sector  of  Figure  6.22 


Figure  6.24:  Active  sector,  zoom  in  of  Figure  6.23 


Parte!  Radiation  Pattern 


Figure  6.25:  Radiation  Pattern  of  Figure  6.22 
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Figure  6.26:  Active  sector  based  on  element  criteria,  amax  =  22°,(#o,  4>o)  =  (37.38°,  36°) 


Panel  Radiation  Pattern 


Figure  6.27:  Array  pattern  pertaining  to  active  sector  of  Figure  6.26  ,  (#0,  <pa)  =  (37.38°,  36°) 


Figure  6.30:  Active  sector  based  on  element  criteria,  amax  =  22°  ,(0q,  <f>0)  =  (47.38°,  36°) 


Panel  Radiation  Pattern 


(180  <  <]>  <  360)  9  (deg)  (0<<|><180) 


Figure  6.31:  Array  pattern  pertaining  to  active  sector  of  Figure  6.30,  (0q,  4>o)  =  (47.38°,  36°) 


Panel  Radiation  Pattern 


Figure  6.35:  Array  pattern  pertaining  to  active  sector  of  Figure  6.34,  ($0,  (po)  =  (57.38°,  36°) 
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Figure  6.37:  Array  pattern  pertaining  to  active  sector  of  Figure  6.36,  (0q,  <t>o)  —  (37.38°,  36°) 
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Figure  6.38:  Active  sector  based  on  panel  criteria,  amax  =  22°, ((90, 0o)  =  (42.38°,  36°) 
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Figure  6.39:  Array  pattern  pertaining  to  active  sector  of  Figure  6.38,  (0q,  4> o)  =  (42.38°,  36°) 


Figure  6.42:  Active  sector  based  on  panel  criteria,  amax  =  22°,(<90,  4>o)  =  (52.38°,  36°) 


Pane!  Radiation  Pattern 


Figure  6.43:  Array  pattern  pertaining  to  active  sector  of  Figure  6.42,  (00, 4>o)  =  (52.38°,  36°) 


Panel  Kadiation  Pattern 


Figure  6.45:  Array  pattern  pertaining  to  active  sector  of  Figure  6.44,  (0q,  4>o)  =  (57.38°,  36°) 
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Appendix  A 

Element  Patterns  in  source  (local) 
coordinates 

A.l  Far-field  of  two  crossed  dipoles 

We  evaluate  the  far-field  due  to  two  crossed  infinitesimal  dipoles,  A/4  above  an  infinite 
conducting  ground.  First  dipole  points  in  the  x-direction  with  current  Iq  while  second  the 
dipole  points  in  y-direction  with  current  Jo  exp(±jn/2).  The  far  field  is 

E(r,  9,  <f>)  =  Eq-  -  sin(^  cos  0)e±J^(0  cos  9  ±  j4>)  (A.l) 

where 

Eo  =  Co^  (A.2) 

and  +/—  denotes,  in  general,  left/right  handed  elliptically  polarized  wave,  and  when  9  =  0, 
LHCP/RHCP  wave. 

In  (A.l),  it  is  convenient  to  replace  sin(7r/2  cos  9)  with  cos^/2)  9.  This  approximation  is 
adequate  for  most  practical  applications  and  it  is  exact  at  9  =  0  and  9  =  7r/3.  Thus, 

E(r,  9 ,  <f>)  =  E0^J^Vcos9e±j<t,(e  cos  9  ±  j$).  (A.3) 

r 

This  expression  can  be  further  simplified  by  assuming  a  CP  wave  for  all  ( 9 ,  <p)  as  follows 

E(r,  9,  <f>)  =  EQ^—Vcos9e±j*(e  ±  j$)  (A.4) 

r 

or  in  the  theoretical  limit  for  isotropic  CP  source 

p—jkr 

E(r,  9 ,  <t>)  =  Eo - e^(d  ±  (A.5) 

r 
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A. 2 


Approximate  expressions  for  array  element  pat¬ 
terns 


For  a  precise,  low  sidelobe  design,  one  needs  to  know  the  radiation  and  impedance  properties 
of  the  array  element  in  a  mutually  coupled  spherical  array  environment.  This  requires  a  rig¬ 
orous  analysis  for  the  specific  element  and  array  geometry  [36,  37].  For  a  sphere  with  large 
radius,  the  local  spherical  array  geometry  can  be  approximated  by  an  equivalent  infinite 
planar  array.  Furthermore,  if  high  accuracy  is  not  desired,  or  for  a  rough  array  perfor¬ 
mance  estimation,  one  can  use  the  single  mode  theory  of  the  uniform  infinite  planar  arrays, 
which  gives  reasonable  but  simple  expressions  for  the  element  magnitude  and  phase  patterns. 
For  a  fundamental  radiator,  two  crossed  dipoles  above  the  ground  plane,  the  approximate 
expression  for  the  far  field  due  to  (n,  m)  element  is 


Enmfo  0)  =  ®n m  0)  (A. 6) 

Trim 

where  anm  are  the  complex  excitation  coefficients  and  e(9,  <j>)  is  the  element  pattern  which 
can  be  written  in  the  following  form: 

e(0, 4>)  =  (1-  I  ra(0o,0o)  |2)  Vcos  9  e±j4>  (6  cos  9  ±  j4>).  (A.7) 

Consequently,  the  respective  element  gain  is 


G(e,4>)  =  G^u  (1-  I  ra(0o,<£o)  |2)  COS0  i(l  +  cos20). 
Each  term  in  the  equations  above  has  a  distinct  physical  meaning: 


GCell 


4  ‘kAcm 
A2 


4A.8). 


(A.9) 


is  the  unit  cell  area  gain,  cos  9  accounts  for  the  reduction  of  the  projected  unit  cell  area  for. 
off-broadside  field  points,  (1—  |  r(f?o,  0O)  |2)  takes'  into  account  the  reduction  of  element  gain 
due  to  active  impedance  input  mismatch  losses  at  scan  angles  (8q,  4> 0)  while  (1  +  cos2  0) /2 
represents  polarization  losses. 
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Appendix  B 

Spherical  Array  Analysis 


B.l  Array  model 

A  typical  spherical  array  is  shown  in  Figure  B.l.  The  antenna  elements  are  placed  on  the 
surface  of  radius  R.  The  distribution  of  the  elements  may  be  general;  however,  in  Figure  B.l 
to  Figure  B.3  we  show  the  arrangement  where  the  elements  are  aligned  along  z=constant 
circles.  The  angular  spacing  in  elevation  between  the  circles  is  uniform.  The  azimuthal 
element  spacing  on  each  circle  is  also  uniform  but  slightly  differ  from  circle  to  circle  due 
to  different  circle  radii.  The  number  of  elements  on  each  circle  was  determined  so  that  the 
azimuth  element  spacing  on  each  circle  is  roughly  the  same  although  periodic.  Each  element 
is  designated  by  two  indeces  (n,  m)  where  n  designates  the  circle  (row)  counting  from  the 
2-axis  and  m  designates  the  element  on  a  nth  circle  counting  from  x-z  plane  according  to 
’’right-hand  screw”  rule  pointing  into  +z  direction.  Therefore  the  (n,m)  =  (0,0)  is  the 
top  element  located  on  2-axis.  Each  element  is  assumed  to  produce  an  elliptically  polarized 
radiation  having  a  maximum  along  its  own  radial  direction.  The  antenna  elements  are  phased 
such  that  the  pattern  produced  by  the  array  has  a  maximum  along  the  ( 90 ,  0O)  direction, 
as  the  position  of  a  typical  element  at  Pnm  on  the  spherical  surface  is  represented  by  the 
coordinates  (R,an,  /3nm)  as  shown  in  Figure  B.4.  The  coordinates  of  a  far-field  point  P  are 
(r,  9,  <p)  and  the  beam  pointing  direction  is  specified  by  (9q,  </>o)  as  illustrated  in  Figure  B.5. 


B.2  Global  and  Local  Coordinate  Systems 

We  define  a  global  rectangular  coordinate  system '(2,  y,  z)  with  unit  vectors  (x,  y,  z),  and  cor¬ 
responding  spherical  coordinates  (r,  9 ,  4>)  with  respective  unit  vectors  (r,  6, 0).  The  far-field 
observation  point  is  P(r,  9,  <j>)  or  in  vector  form  P(r)  where  r  =  r r  as  shown  in  Figure  B.4. 

At  each  element  (n,  m)  we  also  define  the  local  coordinate  system  (xnm,  ynm,  znm)  (drawn 
in  red  color)  with  unit  vectors  (xnm,ynm,znm)  and  corresponding  spherical  coordinates 
(rnmi  9nm,  (pnm)  with  unit  vectors  ( fnm ,  6nm,  4>nm).  In  local  coordinates  the  far-field  observa¬ 
tion  point  is  P(rnm,  9nm,  </>nm),  see  Figure  B.4. 
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B.3  Coordinate  Transformation  Relations 


The  field  produced  by  a  typical  element  at  Pnm,  is  initially  expressed  in  the  local  coordi¬ 
nate  system  (xnm,ynm,  znm),  with  the  origin  at  Pnm  and  the  znm  axis  being  normal  to  the 
spherical  surface  as  shown  in  Figure  B.4.  The  transformation  relations  for  going  from  the 
(nm)th  local  to  global  coordinates  can  be  obtained  by  using  the  general  method  of  orthogonal 
transformation  of  coordinates.  In  general,  (9nm,  4>nm)  can  be  expressed  as  [38]  to  [41] 

COS  9nm  =  f  •  Znm  (B.l) 

cot  4>nm  =  ^  (B.2) 

^  *  Ynm 

where  (xnm,ynm,znm)  are  unit  vectors  in  {xnm,ynm,znm)  directions,  and 

r  =  xsin5cos^  +  ysin5sin0  +  zcos0.  (B.3) 

The  respective  local  unit  vectors  (xnm,  ynm,  znm)  are  related  to  the  unit  vectors  of  the 
global  coordinate  system  (x,  y ,  z)  by  the  following  matrix  equation: 

(B.4) 

L*J 

where  Anm  is  a  rotation  matrix  and  is  given  by 

"  cos  an  cos  (3nm  cos  an  sin  0nm  -  sin  aT 
Anm  =  sin  f3nm  COS  Pnm  0 

_  sin  an  cos  /3nm  sin  an  sin  /?nm  cosan 

After  using  (B.l)  to  (B.5),  the  following  two  relations  are  obtained: 

cos  9nm  =  sin  an  sin  0cos(0  —  0nm)  +  cos  an  cos  9 
cos  an  sin  9  cos(<p  -  Pnm)  —  sin  an  cos  9 

COt  <Pnm  —  t  ^  "  77  7  \  * 

sm  9  sm(<p  -  pnm) 

Similarly,  the  local  unit  vectors  ( 0nm ,  d>nm)  can  be  transformed  into  global  coordinates  via: 


X 

Ynm 

—  1 

—  nnm 

y 

m  ^nm 

i 

(B.6) 

(B.7) 


'  Km  ' 

&nm 

bum  * 

'  6  ' 

.  *Pnm  . 

bnm 

&nm 

.4>. 

where 


^nm. 


_  cos  9  sin  qn  cos(</>  —  /3nm)  —  sin  9  cos  an 


bnm  — 


sin  9nm 

sin  On  sin(^  —  /?nm) 
Sin  9nm 


and 


sin  9nm  =  yfl  —  cos20nm. 

Taking  inverse  of  (B.8)  where  determinant  A  =  a|m  +  =  1?  we  obtain 

&nm  t 
brim  0 


'  8 

®nm 

.  4*nm  . 

(B.8) 

(B.9) 

(B.10) 

(B.ll) 

(B.12) 
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B.4  Element  patterns  in  local  coordinates 


The  far  electric  field  produced  by  a  single  element  at  the  field  point  P(r,  0,  0)  can  be  written 
as  follows 

krn7a) 

Enm(  Tnmi  @nrm  — ®n7n  ($nmj  $nm)  (B.13) 


with 


nm)  —  ®nm^nm(^nm :  $nm  )  +  0  nm^'4>Tim,{@nm'i  4*nm) 


(B-14) 


where  egnm  and  e^,nm  are  the  complex  element  pattern  factors  of  the  two  field  components. 
Basically  they  represent  the  field  produced  by  a  single  element  in  a  match-terminated  spher¬ 
ical  array  environment.  In  general  egnm  and  e^m  are  of  the  form 


eenm(0nm,<t>nm)  =|  e6nm(0nm,  <pnm)  | 


(B.15) 


e^nm(9nmAnm)  =|  e^m(6U,  0™)  |  ej5 .  (B.16) 

Here,  5  is  the  phase  difference  between  the  two  components  of  the  field.  In  general  the  field 
is  elliptically  polarized.  For  a  circularly  polarized  field  e8nm  =  e^m  and  5  =  ±7t/2  where  a 
plus  sign  corresponds  to  RHCP  and  a  minus  sign  to  LHCP  outgoing  wave. 


B.5  Element  patterns  in  global  coordinates 

The  field  in  (B.13)  can  be  expressed  in  global  coordinates  using  the  coordinate  transformation 
relation  (B.8)  as  follows: 


g j(ut-kr) 


Enm(r ,  0,  <f>)  i)  —  enm(0,  (j> ) 

r 

(B.17) 

where 

enm(0, 4>)  =  [0eflnm(0,  0)  +  0e^m(#,  0)]e^"m 

(B.18) 

and  (from  (B. 12))  —  —  -  -- 

0)  —  &nm  ^Onm^nmi  “1"  ^ nm  &4> nmifinm j  finni) 

(B.19) 

4*)  =  bnm  0nm)  flnm  <ftnm) 

(B.20) 

with 

Ip nm  =  kpnm  •  r  =  kR  [sin  an  sin  9 cos (6  -  /3nm)  +  cos  an  cos  9\. 

(B.21) 

To  facilitate  the  computation  of  the  element  excitation  coefficients,  we  rewrite  (B.19)  and 

(B.20)  in 

the  form 

egnm(0,4>)  =|  eenm{0,(l>)  \  e?*™™ 

(B.22) 

e0nm(<9, 0)  =1  e0nrri(0, 0)  |  eWnm(M)  eJl5 

(B.23) 

and  then 

®nm 

(0, 0)  =  [0  |  egnm(d,  0)  |  +  0  |  e*nm(0, 0)  |  e*]  e**». 

(B.24) 
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To  form  a  beam  in  the  (9q,  <Po)  direction,  all  element  contributions  must  coherently  add  up 
in  the  far  zone.  This  means  that  the  phase  of  the  element  excitation  coefficients  must  be 
the  negative  of  the  phase  of  the  enm.  For  CP,  we  must  have  |  egnm(90,  <po)  |=|  e^nm(9o,  <t> o)  |, 
finm(6o,  0o )  =  ¥Wi(#o,  <Po)  and  5  =  ±7t/2.  This  would  require  separate  magnitude  and 
phase/time  controls  for  each  field  component  which  is  not  practical.  A  single  magnitude 
and  phase/time  control  per  element  is  instead  used  resulting  in  some  acceptable  polarization 
loss. 

In  the  above  expressions  egnm  and  e^m  must  be  known  whether  by  EM  simulation,  by 
measurements  or  by  some  kind  of  realistic  approximation. 


B.6  An  Example 

Two  examples  will  be  presented.  In  the  first  example  we  consider  a  CP  isotropic  (nm)th 
array  element  on  the  sphere  with  element  pattern  given  in  local  coordinates  as 

enm(0„m,  <j>nm)  =  E0(6nm  ±  j<pnm)  e±j4>nm.  (B.25) 

In  this  case  we  can  write 

m(9rtm>  *Pnm)  =  1  (B.26) 

C^nm(^nmi  *Pnm)  =  (B.2T) 

Applying  coordinate  transformation  (B.19)  and  (B.20)  yields  - -  _ - i . 

ednm(Q,  <t>)  =  -(flnm  ±  jbnm )  (B.28) 


€<fmm(9.  (flnm  i 

Next,  we  rewrite  (B.28)  and  (B.29)  in  the  form 

eenm(0,<P)  =  -  I  eenm(6,<t>)  \ 

=  ~(±j)  |  e^nje^)  | 


Noting  that 


i  jbnm  j —  1 


(B.31) 


(B.32) 


tan  vnm(6, 4>)  =  tan $nm(0,  <j>)  =  tan  ipnm(9,  <f>)  =  — 

0>nm 

which  after  substituting  (B.9)  and  (B.10)  into  (B.33)  becomes 


tan  vnm(6 ,  <j> ) 


_ sin  an  -  (5nm) _ 

cos  9  sin  an  cos  (<p  —  /3nm)  —  sin  9  cos  an ' 


(B.30)  and  (B.31)  can  be  written  in  the  form 


(B.33) 


(B.34) 


e6nm(9,<P)  =  ~e±ju-^ 


(B.35) 
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r 


) 


* 


•s 


and 

«*.»(«.  i)  =  -(±>)  (B.36) 

The  element  pattern  (B.25)  can  be  now  rewritten  in  global  coordinates  as 

enm{9,  <j>)  =  -E0(6  ±  j$)  e^M)e^nme±j4.nm  (B.37) 

As  a  second  example  we  consider  a  typical  and  simple  approximation  for  eg  and  as 
given  by  (A.3).  Although  the  phase  S  =  ±ir/2  the  wave  is  CP  only  on  the  element  radial 
axis,  i.e.,  when  9nm  =  0.  Thus  we  write 

enm(9nm,  4>nm)  =  Eq  (enm  cos  8nm  ±  j 4>nm)  e±j*nm  (B.38) 


^rwn(^nm)  $nm)  —  COS  9nm 
£(fmm  (^n77i7  4*nm  )  =  ±j 

where  eSnm  =  e^m= 0  for  9nm  >  7r/2. 

Following  the  same  procedure  of  first  example,  we  get 


eenm{BA)  =  -  I  eenm(B,4>)  \ 

etnm(e,  <f>)  =  ~(±j)  |  e*nm(0,  <l>)  I 

where 

I  eenm(Q,<t>)  1=1  cos  9nm  i  jbnm  | 
|  i<j>)  |=|  Q-nm  i  j bnm  COS  9nm  | 


and 


tan  I9nm  — 


tan</?„ 


nm  COS  9nm 
nm  COS  8nm 


Then,  (B.38)  can  be  rewritten  in  global  coordinates  as  . 


(B.39) 

(B.40) 


(B.41) 

(B.42) 

(B.43) 

(B.44) 

(B.45) 

(B.46) 


,(M)  =  -E„  |0  I  etnje,i)  I  ±jj,  I  I 


(B.47) 
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'igure  B.2;  Spherical  array,  side  view 


9 

B.4:  Global  and  local  coordinate  systems 
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Appendix  C 

Maximum  Gain  Theorem  for 
Conformal  Arrays 


The  optimization  problem  for  antenna  arrays  has  been  treated  by  many  authors  in  the  past. 
The  first  solution  was  given  by  [15]  where  it  is  shown  that  in  the  absence  of  mutual  coupling, 
the  maximum  directivity  is  equal  to  the  sum  of  the  directivities  of  the  individual  radiators. 
For  example,  for  arrays  of  N  isotropic  sources  with  element  spacing  d  =  A/2  the  maximum 
directivity  Dmax  =  N.  If  there  is  a  mutual  power  exchange  between  the  array  elements,  as 
shown  in  [15]  the  highest  Dmax  may  be  obtained,  usually  for  element  spacing  d  <  A/2,  and 
in  the  limiting  case  as  d  — >  0,  Dmax  — ►  N2.  These,  so  called  super  gain-arrays,  as  is  well 
known,  require  large  element  currents  distributions,  with  inter-element  phase  close  to  180 
deg  yielding  low  efficiency  and  extremely  narrow  frequency  bandwidth,  and  therefore,  they 
are  not  practical.  Over  the  years,  this  work  has  been  followed  by  many  papers  dealing  with 
the  optimization  of  the  performance  of  arbitrary  antenna  arrays  most  of  them  focusing  on 
end-fire  arrays  [16].  A  notable  theorem  concerning  the  maximum  realized  gain  of  arbitrary 
broadside  arrays  is  proven  in  [17].  The  theorem  basically  states  that  the  maximum  gain  of 
a  nonuniform  conformal  array  is  equal  to  a  sum  of  the  array  element  gain  pattern  values 
in  the  main  beam  direction.  This  is  obtained  with  the  following  characteristic  excitation: 
complex  voltage  excitation  coefficients  are  proportional  to  the  complex  element  (field)  pattern 
values  in  the  main  beam  direction  and  a  phase  factor  accounting  for  the  different  element 
propagation  paths  to  a  common  phase'  front.  The  theorem  does  not  apply  to  end-fire  and 
super-gain  arrays. 

For  convenience,  below  we  re-derive  the  excitation  criteria  for  maximum  gain  of  the 
non-uniform,  an  arbitrary  curved-surface  arrays  originally  presented  in  [17]. 

C.l  The  Theorem 

We  consider  an  arbitrary  array  of  N  elements  as  shown  in  Figure  C.l.  The  array  elements  are 
fed  by  voltage  generators  Vgi,  with  impedance  Zgi,  where  i  =  1, . . . ,  N  denotes  the  element 
serial  number.  The  a*  and  bu  in  Figure  C.l,  are  the  voltages  of  incident  and  reflected  signals 
at  the  array  element  terminals.  The  a*  are  also  called  the  array  excitation  coefficients.  Each 
array  element  is  characterized  by  its  element  gain  pattern  gi(9,(f>),  where  the  angles  0,4> 
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define  the  far  field  observation  direction  in  spherical  coordinates.  It  is  important  to  mention 
that  gi  represents  the  element  gain  pattern  of  the  ith  element  in  a  mutually  coupled  array 
environment  (which  can  be  measured  when  element  i  is  excited  while  all  other  array  elements 
are  match-loaded).  For  simplicity,  it  is  assumed  that  the  array  elements  are  lossless. 

The  Maximum  Gain  Theorem  for  conformal  arrays  states  that  the  maximum  possible 
gain  (for  a  particular  polarization)  of  a  conformal  array  in  the  main  beam  direction  (9q,  <fs0) 
is  the  sum  of  individual  element  gain  values  in  this  direction,  i,e,, 

N 

G(&0,  <Po)maz  =  53  900,  #o)-  (C.l) 

i=l 

This  can  be  achieved  if  the  magnitudes  of  the  excitation  coefficients  are  proportional  to  the 
element  gain  values  in  the  (9q,  <po)  direction,  i.e., 

I  <k  |oc  \fg0o^k)  (C.2) 

and  the  phase  of  the  excitation  coefficients  is  such  as  to  constructively  add  all  element  field 
contributions  in  the  same  direction. 


C.2  Proof 

To  prove  (C.2),  we  consider  an  arbitrary  array  of  N  elements  as  shown  in  Figure  C.l.  In  this 
figure,  r'  is  the  position  vector  of  the  i-th  element  with  respect  to  the  origin  of  the  spherical 
coordinate  system  (r,  9,  <j>).  In  this  case  the  far  field  due  to  element  i  is 

e-ifen 

E0i,  9, 4>)  =  a{e0,  <p) -  (C.3) 

n 

where  a*  axe  incident  complex  excitation  coefficients 

a i  =|  (k  I  ejai  .  (C.4) 


and 

ei(O,<l>y=\e0,<f>)  (C.5) 

are  the  complex  element  patterns.  In  general  the  vector  element  patterns  are  different  for 
each  element  and  represent  the  radiation  pattern  of  a  single  element  in  a  mutually  coupled 
array  environment.  Furthermore  in  (C.3),  is  the  distance  between  the  i-th  element  and 
the  far-field  observation  point  (r,9,(f>).  Consequently,  the  element  gain  is 

gt{9,  <f> )  =  (C.6) 

Co-Poi 


where 


Poi  = 


(C.7) 
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r 


and  Co  =  376.7  Cl  is  the  free  space  impedance.  For  simplicity,  we  assume  that  Zgi  =  l,  thus 
in  (4b),  |  Oi  |2  represents  the  total  incident  (available)  power  in  the  ith  element.  Substituting 
(3a)  into  (4a),  the  element  gain  becomes 

A'it 

a.(M)  -  «(M)  I*  •  (C.8) 

Using  (C.3),  by  superposition,  the  total  array  far  field  is 

p—jkr  N 

E( r)  = - £  a%ex(6,  (C.9) 

r  i=i 


where  f  =  r/r  is  the  unit  vector  toward  the  observation  point  (r,  9,(j>).  In  this  expression  we 
have  used  the  approximation  r{  ~  r  in  the  amplitude  factor  1  /n.  Consequently,  the  array 


gam  is 


where 


G(e,<f>)  = 


E(r)E*(r) 


Po  = 


Co  Po 
1  q«  l2 

Airr2 


(C.10) 

(C.U) 


Substituting  (C.9)  into  (C.10),  the  array  gain  becomes 


G(M) 


47T 1  Ejli  ^(fl,  <ft)eJ'fcr<'ro  |2 
Co  Eili  I  ^  P 


(C.12) 


To  prove  the  maximum  gain  theorem  we  address  the  following  question:  for  a  given  ej, 
what  cij’s  are  required  to  achieve  a  maximum  array  gain  in  the  particular  (00, 4>o)  direction. 
Thus,  we  must  maximize  G(9o,<Po)  as  given  by  (C.12)  with  respect  to  |  a*  |  and  a*.  This  is 
done  in  two  steps. 


step  1:  For  the  moment  we  assume  that  the  |  a,  |  are  known.  In  this  case  in  view  of  (C.4) 
and  (C.5)  we  see  that  (C.12)  is  a  maximum  when 

kv'i  ■  r0 '+  ipi(8o,  <j>o )  +  otl  =  0  (C.13) 

from  where 

OLi  =  —kx[  ■  r0  -  ipi(9 o,  (j)o).  (C-14) 

As  expected,  the  first  term  in  (C.14)  accounts  for  the  different  propagation  path  lengths 
while  the  second  term  takes  into  account  the  different  element  phase  pattern  values  in  the 
(6>0,  <po)  direction.  With  such  a,,  all  array  element  far-field  contributions  will  add  up  in  phase. 


step  2:  With  a,  of  (C.14),  relation  (C.12)  can  be  written  in  the  following  form: 


G(80,<f>0) 


4tt[E^i  1  at  II  ei(6*o,^>o)  ]]2 
Co  E£Li  I  a*  I2 


(C.15) 
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In  order  to  determine  the  j  a*  |  that  gives  the  maximum  G(6 0, <f>o)  it  is  convenient  to  renor¬ 
malize  |  a i  |  according  to 

*  *  =.  (C.16) 


ck  1= 


1  I  °i 


In  this  case  (C.15)  becomes 


4tt  * 

C(«o,*)  =  rE|ai||ei(»o,Ai)  I]2. 

SO  i=i 

By  introducing  the  excitation  vector 

a  =  (|  Si  |,|  a2 d^r  [) 


(C.17) 


(C.18) 


and  the  element  pattern  vector 

©  =  (1  ei($o,0o)  |) |  62(^01^0)  !>•••)!  ^iv(^O)^o)  |)  ..  (C.19) 

(C.17)  can  be  further  rewritten  as  the  complex  scalar  product 

G(90,<l>0)  =  ^[a-e}2.  (C.20) 

so 

Using  Schwartz’s  inequality  |  (a,  e)  |<|  a  |J  e  J,  we  see  that  the  excitation  vector  which  yields 
the  maximum  gain  is  determined  by  |  a  |=  C  |  e  |,  where  C  is  the  constant  to  be  determined. 
In  component  form  this  can  be  written  as 

lai|=C']ei(Mo)|  (C.21) 

where  C  is  the  proportionality  constant  to  be  determined.  Thus, 

4tt  n 

G(i 9o,  (po)max  =  C2-E  I  et(0o,  <f> 0)  |2]2-  (C.22) 

So  i=i 


To  determine  C  we  write 
from  which 


I  |2f=  C2  |  ei(^0,  <f>o)  I2 
I  a i  |2=  C2  J2  !  ei{6 0, 0o)  I2  • 

1=1  i= 1 


With  the  help  of  (C.16)  we  find  that  the  left  hand  side  of  (C.24)  is 

£  I  *  1;-  D-dry  ■  f  -  £  - 1- 

<= 1  »= 1  vE)=i  I  I2  *=1  ^i=i  l  I 

Substitution  of  (C.25)  into  (2.23)  yields 

N 


l  =  C2'£\ei(e0l<t>0) 

t=l 


(C.23) 

(C.24) 

(C.25) 


(C.26) 


s 

# 


S 
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so  that 


C2  = 


£,=i  |  e{(&0, 4>o)  |2‘ 


Therefore,  substituting  (C.27)  into  (C.22),  the  maximum  gain  is 

A-jr  * 

Gfo,  I" 

M  i=l 


or  in  terms  of  the  element  gain  pattern  as  given  by  (C.8) 

N 

G(9o,  <£o)  max  —  !>(*.*>)■ 
i= 1 


(C.27) 


(C.28) 


(C.29) 


As  shown,  this  can  be  achieved  if  the  amplitudes  |  a,  |  and  phases  a;  of  the  complex  excitation 
coefficients  a*  satisfy  (C.21)  and  (C.14),  respectively.  To  elaborate  on  the  former  we  write 


\ai\=C\ei(90,(f>0)  |= 


I  et(flcb  fo)  | 

I  «,(«>,*>)  I2' 


(C.30) 


For  convenience,  we  normalize  (C.30)  with  respect  to  the  first  array  element,  which  yields 


1  %  I  _  I  <k  1  _  1  ft(0 o,0q)  1  ft(0o,<fo) 

I  °i  I  I  I  I  ei(^o,^o)  I  ^  9i(0o,(/)o) 


(CM) 


Thus,  to  obtain  the  maximum  gain  of  the  conformal  array  we  have  found  that  the  magnitude 
of  the  excitation  coefficients  must  be 


Q.E.D. 


I  Oi  1=1  “i  I 


ft  (flo,  0o) 
\  9i(0o,<l>o) 


(C.32) 
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